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Abstract: 

QCD predicts that strongly interacting matter will undergo a transition from a state 
of hadronic constituents to a plasma of unbound quarks and gluons. We first survey the 
conceptual features of this transition and its description in finite temperature lattice QCD, 
before we address its experimental investigation through high energy nucleus-nucleus col- 
lisions. After considering the conditions achievable in such collisions, we discuss the pos- 
sible probes to check if the produced medium in its early stages was indeed deconfined. 
We then elaborate the method that has emerged and the results which were obtained 
using the most extensively studied deconfinement probe, the suppression of charmonium 
production. In closing, we discuss possible supporting information provided through the 
study of soft hadronic probes. 



Contents 



1. Introduction 

2. States of Matter in QCD 

2.1 Prom Hadronic Matter to Quark-Gluon Plasma 

2.2 Critical Behaviour in QCD 

2.3 Dcconfinement and Chiral Symmetry Restoration 

2.4 Deconfinement and Percolation 

3. Conditions in Nuclear Collisions 

3.1 High Energy Collisions and Nuclear Transparency 

3.2 Collision Evolution and Initial Conditions 

3.3 The Onset of Deconfinement 

4. Probes of Primordial Matter 

4.1 Evolution Stages 

4.2 Quarkonium Dissociation 

4.3 Jet Quenching 

5. J/V' Suppression in Nuclear Collisions 

5.1 The Hadroproduction of Charmonium 

5.2 Pre-Resonance Suppression 

5.3 Anomalous Suppression 

5.4 Pt Dependence 

6. Soft Hadronic Probes 

6.1 Retrospective Probes 

6.2 Hadronisation and Freeze-Out 

6.3 Parton Cascade Models 

7. Summary 



1. INTRODUCTION 

Over the past hundred years, our ideas about the ultimate constituents of matter have 
undergone a considerable evolution. Atoms were found to be divisible into electrons and 
nuclei. Nuclei in turn consist of nucleons, bound together by strong short-range forces. 
With the advent of the basic theory of strong interactions, quantum chromodynamics 
(QCD), has come the conviction that nucleons - and more generally, all strongly inter- 
acting elementary particles (hadrons) - are bound states of quarks. Quarks are point-like 
and confined to "their" hadron by a binding potential Vo(r) which increases linearly with 
quark separation r, 

^o(r) ~ ar, (1) 

where the string tension a measures the energy per unit separation distance. Hence an 
infinite amount of energy would be needed to isolate a quark; it cannot exist by itself, and 
it is therefore not possible to split an isolated hadron into its quark constituents. The 
elementary particles of strong interaction physics remain elementary in the sense that 
they are the smallest entities with an independent existence; however, they have become 
composite, since they are bound states of quarks. What effect will this have on the states 
of strongly interacting matter? 

To get a first idea of what the quark infrastructure of elementary particles implies for 
the behaviour of matter at extreme density, consider a very simple picture. If nucleons, 
with their intrinsic spatial extension, were both elementary and incompressible, then a 
state of close packing would constitute the high density limit of matter (Fig. |l]a). On 
the other hand, composite nucleons made up of point-like quarks will start to overlap 
with increasing density, until eventually each quark finds within its immediate vicinity a 
considerable number of other quarks (Fig. |I]b). It has no way to identify which of these 
had been its partners in a specific nucleon at some previous state of lower density. Beyond 
a certain point, the concept of a hadron thus loses its meaning; at extreme density, we 
are quite naturally led to a medium whose basic constituents are unbound quarks. 




a b 

Figure 1: Strongly interacting matter as nuclear matter at a density of closely packed 
nucleons (a) and as quark matter at much higher density (b). 

Ever since the big bang, our universe has been expanding; if we could reverse time 
and let the film run backwards, we would find matter of higher and higher densities. Up 
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to nuclear density, we have experimental information to guide our understanding of such 
matter. What happens beyond this point? This question, the nature of matter at extreme 
density, will be the central topic of this report. The primordial matter in the very early 
universe provides a particularly striking instance of such a system, but not the only one, 
as we shall see shortly. 

On the theoretical side, studies of statistical QCD, in particular through computer 
simulations, have confirmed the transition from hadronic matter to a plasma of unbound 
quarks and gluons. This led to the challenge to produce and study the predicted new 
state of matter in the laboratory. The only way to achieve this, as far as we know, is to 
collide two heavy nuclei and study the resulting small and short-lived droplets of hot and 
dense medium. For such studies it is essential to have viable probes which can determine 
if the droplets in the early stages of nuclear collisions indeed consisted of unbound quarks 
and gluons: we have to find signatures of colour deconfinement. 

In the next section, we shall survey concepts and theoretical results concerning the 
states of matter in strong interaction physics. Following that, we shall address the ther- 
modynamic conditions which can be achieved in nuclear collisions. This will set the stage 
for the main problem: probing experimentally whether high energy nuclear collisions have 
produced deconfined matter and studying the properties of the new phase. Finally we 
shall consider the possible support that soft hadronic observables from later stages can 
provide to corroborate the conclusions obtained using early hard probes. 

In view of the many-faceted nature of high energy nuclear physics, it is perhaps ap- 
propriate to emphasize here that the aim of this survey is colour deconfinement and its 
manifestation in such collisions. Many other interesting aspects of the field will unfortu- 
nately but unavoidably have to suffer from this restriction. 



2. STATES OF MATTER IN QCD 

2.1 From Hadronic Matter to Quark-GIuon Plasma 

Confinement is a long-range feature, which prevents the isolation of a single quark. In 
the high density situation of Fig. |l]a each quark finds very close to it many others and 
is thus far from isolated. To see how the short-range nature of dense matter overcomes 
confinement, recall the effect of a dense medium on electric forces. In a vacuum, two 
electric charges eg interact through the Coulomb potential 



where r again denotes the separation distance. In a dense environment of many other 
charges, the potential becomes screened. 




(2) 




(3) 
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where r^) = is the Debye screening radius of the medium; it decreases as the charge 
density of the medium increases. Thus the potential between two test charges a fixed 
distance apart becomes weaker with increasing density. This occurs because the other 
charges in the medium partially neutralize the test charges and thereby shorten the range 
of the interaction. If a bound state, such as a hydrogen atom, is put into such a medium, 
the screening radius will for sufficiently high density become less than the binding 
radius of the atom. Once r£, <^ r^, the effective force between proton and electron 
has become so short-ranged that the two can no longer bind. Thus insulating matter, 
consisting of bound electric charges, will at sufficiently high density become conducting: 
it will undergo a phase transition [|l|, in which charge screening dissolves the binding of 
the constituents, leading to a plasma of unbound charges as a new state of matter. 

The interaction of quarks in QCD is based on their intrinsic colour charge, and in a 
dense medium this charge can be screened in much the same way as an electric charge. 
Hadrons are colour-neutral bound states of coloured quarks; hence dilute hadronic matter 
is a colour insulator. At sufficiently high density, however, we expect colour screening to 
set in, so that the potential (|ip becomes [| 

1 -exp(-^r) l _ 
fir 

as above, the colour screening mass fj, here is also the inverse of the screening radius 
for colour charges. The resulting exponential damping of the binding force will remove 
all long range effects and in a sufficiently dense medium 'melt' hadrons (see Fig. |^) just 
as Debye screening dissociated hydrogen atoms. Colour screening will thus transform a 
colour insulator into a colour conductor, turning hadronic matter into a quark plasma |^. 
The transition from insulator to conductor by charge screening is a collective effect, so 
that we expect a phase transition at the point of plasma formation. 
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Figure 2: Colour screening of the confining potential. 

^The functional form of the screening depends on the form of the unscreened potential this leads 
to the difference between Eqs. (^) and 



V{r) ~ ar 
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When atomic matter is transformed from an insulator into a conductor, the effective 
mass of the conduction electrons is changed. In insulators, the electrons occur with their 
physical mass; in the conducting phase, however, they acquire a different, effective mass, 
due to the presence of the other conducting electrons, the periodic field of the charged 
ions and the lattice vibrations. All these effects combine to produce a mean background 
field quite different from the vacuum. As a result, the insulator-conductor transition 
is accompanied by a mass shift for the electrons. Similarly, the effective quark mass is 
expected to change between the confined and the deconfined phase. When confined in 
hadrons, the basic quarks 'dress' themselves with gluons to acquire an effective constituent 
quark mass of about 300 MeV (1/3 of the proton or 1/2 of the p-meson mass). On the 
other hand, the basic bare quarks in the QCD Lagrangian are almost massless, so that the 
mass of the constituent quarks in the confined phase must be generated spontaneously 
through the confinement interaction. Hence it is likely that when deconfinement occurs, 
this additional mass is 'lost' and the quarks revert to their intrinsic 'bare' mass. 

A Lagrangian with massless fermions - the limiting case of the light up and down 
quarks in the physical Lagrangian - possesses chiral symmetry; this allows a decomposition 
of the quarks into independent left- and right-handed massless spin one-half states, which 
for massive fermions become mixed. For massless quarks, confinement must thus lead 
to spontaneous breaking of chiral symmetry, deconfinement to its restoration. Hence the 
mass shift transition in QCD is often referred to as chiral symmetry restoration. It can, 
but does not need to coincide with deconfinement: when the hadrons are dissolved into 
quark constituents, the liberated and hence now coloured quarks may still interact to form 
coloured bound states. Thus at low temperature and high density, the quark triplets in 
nucleons, once deconfined, might choose to recombine into massive coloured quark pairs 
("diquarks") 0], similar to Cooper pairs in QED. When the density is increased further, 
the diquarks would break up into the massless basic quarks. This results in a three- 
phase picture of strongly interacting matter, with hadronic matter as confined phase, 
then deconfinement, followed by a phase consisting of massive coloured diquark systems, 
and finally, after chiral symmetry restoration, a plasma of coloured massless quarks and 
gluons. Such a three-phase structure would correspond to the succession of insulator, 
superconductor and conductor in atomic matter. 

In relativistic thermodynamics, higher densities can be obtained either by increasing 
the net baryon number density, or by 'heating' the system, so that collisions between its 
constituents produce further hadrons. This leads to the phase diagram shown in Fig. |^: 
for low values of the temperature T and the baryochemical potential hb, associated to the 
net baryon number density, we have confinement and hence hadronic matter; for high T 
and/or fiB, deconfinement sets in and we get a quark-gluon plasma. In the compression 
of cold nuclear matter, by increasing fiB for T = 0, an intermediate diquark phase could 
occur, similar to the superconducting phase in atomic matter. By increasing T at /x^ = 0, 
we are heating mesonic matter until it becomes a quark-gluon plasma. Strong interation 
thermodynamics thus predicts the existence of new, deconfined states of matter at high 
temperatures and densities. 

Before turning to the theoretical basis of these conceptual considerations, we illustrate 
the transition from hadronic matter to quark-gluon plasma by a very simple model. For 
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Figure 3: The phase diagram of strongly interacting matter. 



an ideal gas of massless pions, the pressure as function of the temperature is given by the 
Stefan-Boltzmann form 
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(5) 



where the factor 3 accounts for the three charge states of the pion. The corresponding 
form for an ideal quark-gluon plasma with two flavours and three colours is 
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(6) 



In Eq. the first temperature term in the curly brackets accounts for the two spin and 
eight colour degrees of freedom of the gluons, the second for the three colour, two flavour, 
two spin and two particle-antiparticle degrees of freedom of the quarks, with 7/8 to obtain 
the correct statistics. The bag pressure B takes into account the difference between the 
physical vacuum and the ground state for quarks and gluons in a medium. 

Since in thermodynamics, a system chooses the state of lowest free energy and hence 
highest pressure, we compare in Fig. ^ the temperature behaviour of Eq's. (^) and (^. 
Our simple model thus leads to a two-phase picture of strongly interacting matter, with 
a hadronic phase up to 

Te= (^) '^'5^/^^0.72 (7) 

and a quark gluon plasma above this critical temperature. From hadron spectroscopy, 
the bag pressure is given by B^/^ ~ 0.2 GeV, so that we obtain 

Tc ~ 145 MeV (8) 

as the deconfinement temperature. In the next section we shall find this simple estimate 
to be remarkably close to the value obtained in lattice QCD. 

The energy densities of the two phases of our model are given by 



TT 
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a b 

Figure 4: Pressure (a) and energy density (b) in a two-phase ideal gas model. 



and 

e,, = 37^T' + B. (10) 

By construction, the transition is first order, and the resulting temperature dependence 
is shown in Fig. ^d. We note that even though both phases consist of massless non- 
interacting constituents, in the quark-gluon plasma 

{e-3P) = AB (11) 

does not vanish, as consequence of the difference between physical vacuum and in-medium 
QCD ground state 0. 

2.2 Critical Behaviour in QCD 

Here we shall derive the conceptual considerations of the last section from strong 
interaction thermodynamics, based on QCD as the input dynamics. QCD is defined by 
the Lagrangian 

^ = -^C^r - Y.^^{^r^, + mJ-grA,r%, (12) 

with 

i^;. = (d.At - d^Al - 9flA\Al) . (13) 

Here denotes the gluon field of colour a (a=l,2,...,8) and ipl the quark field of colour 
a (a=l,2,3) and flavour /; the input ('bare') quark masses are given by nif. With 
the dynamics thus determined, the corresponding thermodynamics is obtained from the 
partition function, which is most suitably expressed as a functional path integral, 

Z{T,V) = J dA d^p d^j exp (- j^dPx j^J^ dr CiA.il),^)^ (14) 



since this form involves directly the Lagrangian density defining the theory. The spatial 
integration in the exponent of Eq. ([I^) is performed over the entire spatial volume V of 
the system; in the thermodynamic limit it becomes infinite. The time component Xq is 
"rotated" to become purely imaginary, r = ixo, thus turning the Minkowski manifold, on 
which the fields A and ip are originally defined, into a Euclidean space. The integration 
over r in Eq. (p!4D runs over a finite slice whose thickness is determined by the temperature 
of the system. 

From Z{T, V), all thermodynamical observables can be calculated in the usual fashion. 
Thus 

^ = iTW) ( ^1 (15) 



dT 



V 



gives the energy density, and 

fd\nZ\ , , 

P = ^[-W-), (16) 

the pressure. For the study of critical behaviour, long range correlations and multi-particle 
interactions are of crucial importance; hence perturbation theory cannot be used. The 
necessary non-perturbative regularisation scheme is provided by the lattice formulation of 
QCD 1^; it leads to a form which can be evaluated numerically by computer simulation 0. 
Unfortunately, this method is so far viable only for the case of vanishing baryon number 
density, so that results are today available only for that part of the phase diagram in Fig. 

i 

Without going into the details of computer simulation of finite temperature QCD (for 
a survey, see [^), we summarize what has been obtained so far. The first variable to 
consider is the deconfinement measure given by the Polyakov loop 10 



L(r)~;im exp{-r(r)/r} (17) 

where Vir) is the potential between a static quark-antiquark pair separated by a distance 
r. In the limit of large input quark mass, V"(oo) = oo in the confined phase, so that 
then L = 0. Colour screening, on the other hand, makes Vir^ finite at large r, so that in 
the deconfined phase, L does not vanish. It thus becomes an 'order parameter' like the 
magnetisation in the Ising model; this is zero at high temperatures, but becomes non- 
zero below the Curie point, when the 'up-down' Z2 symmetry of the Ising Hamiltonian is 
spontaneously broken. In the large quark mass limit, QCD reduces to pure S'f/(3) gauge 
theory, which is invariant under a global ^3 symmetry. The Polyakov loop provides a 
measure of the state of the system under this symmetry: it vanishes for Z3 symmetric 
states and becomes finite when Z3 is spontaneously broken. Hence the critical behaviour 
of the Ising model (and more generally, spin systems) as well as that of SU{N) gauge 
theories is based on the spontaneous symmetry breaking of a global Zj^j symmetry . 

For finite quark mass rriq, V{r) remains finite for r —>■ 00, since the 'string' between 
the two colour charges 'breaks' when the corresponding potential energy becomes equal 
to the mass of the lowest hadron; beyond this point, it becomes energetically more 
favourable to produce an additional hadron. Hence now L no longer vanishes in the 
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confined phase, but only becomes exponentially small there, 



L{T) ~ exp{-M,/T}; (18) 

here is of the order of the p-mass, so that L ~ 10~^, rather than zero. Decon- 
finement is thus indeed much like the insulator-conductor transition, for which the or- 
der parameter, the conductivity <j(T), also does not really vanish for T > 0, but with 
cr(T) ~ exp{—AE/T} is only exponentially small, since thermal ionisation (with ioni- 
sation energy AE) produces a small number of unbound electrons even in the insulator 
phase. 

Fig. ^ shows recent lattice results for L{T) and the corresponding susceptibility Xl{T) ~ 
(L^) — (L)^ The calculations were performed for the case of two flavours of light 



quarks, using a current quark mass about four times larger than that needed for the phys- 



ical pion mass [0. We note that L{T) undergoes the expected sudden increase from a 
small confinement to a much larger deconfinement value. The sharp peak of Xl{T) defines 
quite well the transition temperature Tc, which we shall shortly specify in physical units. 

The next quantity to consider is the effective quark mass; it is measured by the ex- 
pectation value of the corresponding term in the Lagrangian, {i/jip) (T) . In the limit of 
vanishing current quark mass, the Lagrangian becomes chirally symmetric and {il)il)){T) 
the corresponding order parameter. In the confined phase, with effective constituent 
quark masses Mq ~ 0.3 GeV, this chiral symmetry is spontaneously broken, while in 
the deconfined phase, at high enough temperature, we expect its restoration. In the real 
world, with finite pion and hence finite current quark mass, this symmetry is also only 
approximate, since {ipilj){T) now never vanishes at finite T. 

The behaviour of {ip'ilj){T) and the corresponding susceptibility Xm ~ d{'ipip) /drUg 
are shown in Fig. |] calculated for the same case as above in Fig. |. We note here 
the expected sudden drop of the effective quark mass and the associated sharp peak in 
the susceptibility. The temperature at which this occurs coincides with that obtained 
through the deconfinement measure. We therefore conclude that at vanishing baryon 
number density, quark deconfinement and the shift from constituent to current quark 
mass coincide. In section 2.3, we shall consider why this is so. 

We thus obtain for = a rather well defined phase structure, consisting of a 
confined phase for T < T^, with L{T) ~ and {iljtp){T) ^ 0, and a deconfined phase for 
T > Tc with L(T) ^ and {il)ip){T) ~ 0. The underlying symmetries associated to the 
critical behaviour at T = Tc, the symmetry of deconfinement and the chiral symmetry 
of the quark mass shift, become exact in the limits — oo and iriq — > 0, respectively. 
In the real world, both symmetries are only approximate; nevertheless, we see from Figs. 
^ and ^ that both associated measures retain an almost critical behaviour. 

Next we turn to the behaviour of energy density e and pressure P at deconfinement. 
In Fig. 1^ it is seen that e/T^ changes quite abruptly at the above determined critical 
temperature T = Tc, increasing from a low hadronic value to nearly that expected for an 



ideal gas of quarks and gluons ||T^. Besides the sudden increase at deconfinement, there 
are two further points to note. In the region Tc<T <2 T^, e and 3P do not yet coincide, as 
also found in the simple model of the previous section; much of the difference is presumably 
due to the difference between physical vacuum and QCD in- medium ground state 0. 



8 



I — ^ 
5 20 



5 , 30 





5 40 



6/g 



Figure 5: The temperature dependence of the Polyakov loop and the associated suscepti- 
bihty in two-flavour QCD pJ]) . 
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Figure 6: The temperature dependence of the chiral condensate and the associated sus- 
ceptibihty in two-flavour QCD . 
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Figure 7: The temperature dependence of energy density and pressure in two-flavour QCD 
fill. 



Furthermore the thermodynamic observables do not quite reach their Stefan-Boltzmann 
values (marked "SB" in Fig. 0) even at very high temperatures. These deviations from 
ideal gas behaviour can be expressed to a large extent in terms of effective 'thermal' 
masses nith of quarks and gluons, with mth ~ gT - |T^. 



Finally we turn to the value of the transition temperature. Since QCD (in the limit of 
massless quarks) does not contain any dimensional parameters, can only be obtained 
in physical units by expressing it in terms of some other known observable which can also 
be calculated on the lattice, such as the p-mass, the proton mass, or the string tension. In 
the continuum limit, all different ways should lead to the same result. Within the present 
accuracy, they define the uncertainty so far still inherent in the lattice evaluation of QCD. 
Using the p-mass to fix the scale leads to Tc ~ 0.15 GeV, while the string tension still 
allows values as large as ~ 0.20 GeV. This means that energy densities of some 1-2 
GeV/fm^ are needed in order to produce a medium of deconfined quarks and gluons. 

In summary, finite temperature lattice QCD at /if, = shows 

• that there is a deconfinement transition coincident with an associated shift in the 
effective quark mass at ~ 0.15 - 0.20 GeV; 

• that the transition is accompanied by a sudden increase in the energy density ( "la- 
tent heat of deconfinement") from a small hadronic value to a much larger value 
near that of an ideal quark-gluon plasma. 

From ab initio QCD calculations we thus have a relatively good understanding of the 
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phase structure expected for strongly interacting matter at vanishing baryon density. 
Two aspects call for some additional discussion. What relation between deconfinement 
and chiral symmetry restoration makes the two phenomena coincide in finite temperature 
QCD? And how does a hadronic medium turn into a quark-gluon plasma on a microscopic 
level - do hadrons somehow fuse to form a deconfined plasma? These questions will be 
addressed in sections 2.3 and 2.4. 



2.3 Deconfinement and Chiral Symmetry Restoration 

The relation between these two transition phenomena has been quite puzzling all 
along; why should they (at least for /is=0) occur at the same temperature? We shall here 



give a somewhat speculative conceptual answer [|18|-||21||. 

In the presence of light quarks of mass rriq, the string between a static quark-antiquark 
pair could in principle break when V{r) = ar = 2mg, leading to 

L{T) ~ (19) 

for the Polyakov loop. For sufficiently small rriq, all variations of L{T) with T would then 
be washed out and the deconfinement transition should disappear ll^. Such an effect 
would be the counterpart of that caused by an external magnetic field H in the Ising 
model. When H = 0, spontaneous magnetisation sets in at the Curie point Tf.; but for 
if 7^0, the field always aligns spins and makes the magnetisation non-zero for all T (Fig. 

We might thus expect that H ~ l/rrig constitutes an effective external field in QCD. 

However, it turns out that this is not the case; the Polyakov loop shows strong varia- 
tions with T even for m^— >0, as seen in Fig. ^ |TT[. Conceptually, there is a good reason 



for this: chiral symmetry breaking in the confined phase keeps the effective constituent 
quark mass Mq finite even for vanishing rriq, and to break the string through formation 
of a normal meson requires 2Mq, not 2mq. The relevant field should thus be 118|-[21 



H ~ 1/Mq, (20) 

where Mq is the constituent quark mass in the relevant temperature region. Hence the 
temperature dependence of Mq plays a crucial role for the temperature dependence of the 
Polyakov loop. 

Putting together what we know about the effective quark mass under different condi- 
tions, we expect the external field 

H ~ l/{mq + Ci(V'V')'/' + T), (21) 

to be determined by a combination of bare quark mass, chiral mass and thermal mass 
(Fig. ITO). For rrig— >oo, we obviously recover the pure gauge case. For small or vanishing 



ruq, {ijj^jj) is the crucial term, since the thermal mass at low T is very small. Below the 
chiral symmetry restoration point T^, the large value of (ipip) keeps the external field H 
small, allowing the Polyakov loop distribution to remain essentially disordered, so that 
the system is in a confined state. At chiral symmetry restoration, (ipip) vanishes and 
hence H suddenly becomes very large; it aligns the Polyakov loops and thereby explicitly 
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Figure 8: Magnetisation in the Ising model with and without external magnetic field H. 
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Figure 9: The temperature dependence of the Polyakov loop for (a) pure gauge theory, 
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Figure 10: The expected temperature dependence of the effective external fleld in full 
QCD. 
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breaks the Z3 confining symmetry. As a result, L increases sharply and deconfinement 
sets in. 

We thus have essentially two limiting situations: 

• For large m^, H is always very small, so that deconfinement is determined by an 
approximate spontaneous symmetry breaking at T = T^. 

• For small m^, H becomes large when {ipip) vanishes; H then aligns the Polyakov 
loops, resulting in deconfinement at T = T^. 

In the second case, we have assumed that chiral symmetry restoration occurs at a tem- 
perature less than the deconfinement temperature in pure gauge theory, as is in fact 
the case, with Tc[^f/(3)] ~ 0.26 GeV and T^{Nf = 2) ~ 0.15 GeV. 

Deconfinement is thus due to spontaneous symmetry breaking only in the large 
rriq limit. For small m^, that symmetry is broken explicitly by an effective external field 
H which becomes large when (ipip) 0. Hence deconfinement and chiral symmetry 
restoration coincide. 

2.4 Deconfinement and Percolation 



Conceptually, the transition from hadronic matter to quark-gluon plasma seems quite 
transparent. As we had seen in the introduction, deconfinement is expected when the 
density of quarks and gluons becomes so high that it no longer makes sense to partition 
them into colour-neutral hadrons, since these would strongly overlap. Instead we have 
clusters much larger than hadrons, within which colour is not confined; deconfinement is 
thus related to cluster formation. This is the central topic of percolation theory, and hence 



a connection between percolation and deconfinement seems very likely ||20| , |22| , |23| . In 
percolation theory, one considers a system of extended geometric objects, such as spheres, 
which can overlap; percolation is said to occur when, with increasing density, a connected 
cluster begins to span the entire system . If the interior of the spheres is coloured, such 
a cluster signals the onset of colour conductivity and hence of deconfinement. Moreover, 



for any value of the quark mass ||2^ , the percolation transition is based on genuine singular 
behaviour as function of the temperature. 

Already the first studies ^2], ^ had noted that percolation of spheres of hadronic size 
cannot mean deconfinement. The percolation density Up for spheres of radius r is given 
by H 

0.34 , , 

(47r/3) r-^ 

Hence for a nucleon radius of 0.8 fm, the critical value is 0.16 fm~'^ and thus just stan- 
dard nuclear density. The percolation of nucleons therefore corresponds to the formation 
of nuclear matter; for deconfinement, the true coloured 'interiors' of the hadrons must 
overlap. These are more easily defined for mesons. 

In QCD, two fundamental triplet colour charges are connected by a string, with 

a ~ 0.16 GeV^ = 0.8 GeV/fm (23) 

as string tension. The transverse radius of these strings is found to be about 0.2 - 0.3 
fm in the range 0.5 < / < 2 fm p6| , p7| , where / is the length of the string. A meson 
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is thus like an ellipsoidof major axis r = 1/2 and minor axis r^, randomly oriented and 
oscillating in time about its center of mass to fill a sphere of radius r. The percolation 
threshold for ellipsoids of aspect ratio r/r^ = 4 is found to be 



0.21 

(47r/3) r rf 



"^p - /._/ox ^ (24) 



for r = 1 fm and = 0.25 fm, this gives 0.8 fm ^ as critical density. With nth = 2ra for 
the 'generic' hadron mass in such a string model, we obtain 

ep ~ 1.3 GeV/fm^ (25) 

as the energy density at the percolation point, in good agreement with the deconfinement 
value (1-2 GeV/fm^) obtained from lattice QCD in section 2.2. Thus deconfinement 
indeed appears to set in at the percolation point for extended hadronic objects; for mesons, 
these are quark- ant iquark strings. For nucleons, the quark triplet structure makes the 
geometry more complex, but conceptually there should be no difference. 

The connection between percolation and deconfinement is presently also under inves- 
tigation on a more formal level, based on the relation between spin and SU(N) gauge 
field systems []TB[ and on the percolation of Polyakov loops [^, As mentioned, the 
introduction of dynamical quarks in QCD destroys the underlying symmetry; the ex- 
pectation value of the Polyakov loop L{T) now never vanishes, so that L{T) can no longer 
be taken as an order parameter to identify the transition. The analogue in the Ising model 
is, as we saw, the introduction of an external field, which also breaks the Z2 symmetry 
and makes the magnetization finite for all temperatures. Although now no thermal phase 
transition can take place, the percolation transition between a region of finite clusters of 
parallel spins to one of infinite clusters continues to exist. Its counterpart in SU(N) gauge 
theory may therefore constitute a general order parameter for deconfinement in full QCD, 
with dynamical quarks [EO . 



3. CONDITIONS IN NUCLEAR COLLISIONS 

3.1 High Energy Collisions and Nuclear Transparency 

Quark matter presumably made up the early universe, and perhaps neutron stars have 
quark cores. But the big bang was long ago, and neutron stars are far away. The rapid 
growth which the field has experienced in the past two decades was to a very large extent 
stimulated by the idea of using high energy nuclear collisions to produce in the laboratory 
droplets of strongly interacting matter large enough and long-lived enough for a study of 
statistical QCD. 

Normal nuclear matter consists of nucleons of mass 0.94 GeV and has a density of 
0.17 nucleons/fm^; hence its energy density is 0.16 GeV/fm^, i.e., well below the 1-2 
GeV/fm^ necessary for deconfinement. We want to collide two heavy nuclei with the aim 
of increasing the density of matter as far as possible beyond this value. To what extent 
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can this be achieved? In a low energy nucleus-nucleus collision, the two nuclei will remain 
intact and simply "bounce off" each other. With increasing energy, they will penetrate 
each other more and more, leading to highly excited nuclear matter, which rapidly breaks 
up into nuclear fragments and some additional mesons. 

If the collision energy is increased still further, nuclear transparency begins to set in: 
the two colliding nuclei pass through each other, leaving behind them a "vapour trail" 
of deposited energy, which eventually decays into hadrons. Different droplets of this trail 
move at different speeds, from slow ones in the center of mass to the fast remnants of 
target and projectile. We refer to elements in the trail as 'droplets', keeping in mind, 
however, that there really is a continuous fluid of such droplets, without any gaps. Since 
nuclear transparency clearly limits the local (droplet) energy density available for the 
medium produced in the collision, its understanding is necessary in order to estimate 
the conditions which can be achieved experimentally. It is a phenomenon which occurs 
already in the collision of two nucleons, which at high energy also pass through each other, 
creating a cascade of secondary hadrons. Most of these are produced outside the primary 
interaction region and long after the primary collision . 

We therefore first ask why in the more elementary process of a nucleon-nucleon collision 
the two collision partners do not stop each other. Nucleons have a spatial extension of 
about 1 fm, and in a nucleon-nucleon collision it therefore takes a certain time tq before 
the entire projectile nucleon has realised that it has hit something. In the rest-frame of 
the projectile, a time tq ~ 1 fm is required to transmit the information about the collision 
from one side of the nucleon to the other. In the rest-frame of the target nucleon, this 
time is dilated to 

to = To (—) , (26) 
Km J 



where Pq = \^P^ + m? is the energy of the incident projectile nucleon, P its lab momen- 
tum and m its mass. Thus at high collision energies, to ^ 1 fm; since the target also has 
a size of about 1 fm, the projectile has long left the target region before it fully realises 
that it was hit. It also takes a time of order 1 fm to form a hadronic secondary of radius 
Rh — 1 fm, like a pion, and hence only very slow secondaries are produced inside the 
target region. For faster ones, the formation time is again dilated, and they are therefore 
formed outside the target region. The production process in nucleon-nucleon collisions, 
beginning about 1 fm after the collision with soft secondaries produced in the target re- 
gion and continuing with further production of faster secondaries later and outside the 
target, is therefore often referred to as inside-outside cascade. 

One thus encounters transparency and a vapour trail of droplets of deposited energy 
already in high energy nucleon-nucleon collisions. Here one moreover observes that the 
nature of all the droplets of different rapidities is more or less the same in their respective 
rest-frames: they decay into the same types of hadrons, and these have the same momen- 
tum distributions, apart from possible effects due to the motion of the droplet In 
other words, the properties of any segment of the vapour trail, measured in its own rest 
system, do not depend on its relative velocity with respect to the others. We will thus 
assume [Q, ^ that the vapour trail is invariant under Lorentz transformations along 
the collision axis. Evidently there are limits to the validity of such a picture: it must 
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break down at both ends of the trail, when the kinematic hmit is reached. Moreover, at 
very high energies, the overall baryon number is negligible in the main (central) part of 
the trail, since the fast target and projectile nucleons emerge as leading particles at high 
rapidity. Hence the droplets have a non-zero baryon number. 

Let us make this more precise. A particle of mass m, moving with a momentum 
p along the beam axis, has a rapidity y = sinh~^(p/m); for the moment, we assume 
vanishing transverse momentum. Since the transformation from one reference system to 
another is just a displacement in rapidity, boost invariance along the beam axis means 
that a physical quantity, such as the deposited energy in the restframe of the droplet, 
cannot depend on the rapidity y. Since the deposited energy appears later in the form 
of hadrons, the number of hadronic secondaries (dNh/dy) produced in a given rapidity 
interval dy must be independent of y, 



'dNh 

. dy 



a(s), (27) 



where a(s) depends on the (squared) collision energy s only. The central rapidity plateau 
defined by Eq. (|27|) is a reasonable approximation to what one observes in high energy 
p — p and p — p collisions. 

We now include transverse momenta and note that the transverse momentum {px) dis- 
tribution of the produced hadrons must also be independent. The full momentum dis- 
tribution of secondaries of species i thus becomes (assuming azimuthal symmetry around 
the collision axis) 

Here ai{s) determines the relative abundance of species i at the given s, and fi{pT,s) 
is the corresponding rapidity-independent transverse momentum distribution, normalized 
to unity. In principle, fi{pT,s) could depend both on the collision energy and on the 
species of hadron involved. One finds, however, that in very good approximation it is 
energy-independent, and the species-dependence enters only through the transverse mass. 



Mpt, s) = f{pT,mi 



2{1 + miX) 



exp{-X [rriTimi) - rrii]}, (29) 



with a universal (i.e., species-independent) constant A^^ ^ 0.15 GeV. Since [mTirrii) —rrii] 
is just the transverse kinetic energy of species i, fi{pT,s) ~ exp{— A Ei} provides a 
universal cut-off in the transverse kinetic energy for all species. This phenomenon is 
generally referred to as m-r-scaling. It is very reminiscent of a Boltzmann distribution 
and has led to the formulation of thermodynamic models for such production processes, 
with A^^ = Th playing the role of a temperature [0. A droplet thus freezes out at 



temperature Th by turning into freely flowing (non-interacting) hadrons. From Eq. 
the average transverse mass of a secondary is 



{mT{mi)) = J 



2(l + miA) 



(30) 
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With A ~ 0.15 GeV, this gives for pions (mT(7r)) ~ 0.37 GeV, for nucleons {mT{N)) ~ 
1.11 GeV. The corresponding average transverse momentum is 0.34 GeV for pions, 0.59 
GeV for nucleons. We thus note that rriT scahng (or thermal production, if we interpret 
it in this way) provides higher average transverse momenta for heavier secondaries. 

The possible thermal nature of hadron production in p — p collisions has recently 



received further support in analyses of species abundances |38|. Experimentally, the 
production rates for up to thirty different species of hadrons and hadron resonances have 
been measured, for different collision energies. The simplest thermal picture has chemical 
freeze-out occurring when a hadron gas, whose constituents are all possible hadron and 



hadron resonance states, reaches a certain temperature ||39|. This would mean that the 
relative abundances of all species at all (high) collision energies are determined by just 
this one parameter. This is in fact not possible for p—p collisions, because of the so-called 
strangeness suppression; the relative abundance of secondaries of non-zero strangeness is 
considerably smaller than predicted. A full description does become possible, however, if 
one further parameter 7^ is introduced to take this into account [|^. It is chosen so that 
7s = 1 implies a fully thermal abundance spectrum, 7^ < 1 strangeness suppression. One 
then finds that for Th — 170 MeV and 7^ ~ 0.5, it is possible to reproduce quite well the 
observed secondary hadron abundances in p — p and p — p collisions from y/s = 20 to 900 



GeV m 



The origin of such thermal behaviour remains rather unclear, all the more so since it 
holds for hadron production in high energy e~^e~ annihilation as well [ETl. It can certainly 



not arise from rescattering among the produced secondaries, since neither the life-time 
nor the density of the produced medium are high enough for this. It may be due to a 
statistical excitation of the vacuum due to the passage of one or more colour charges, 
occurring already on the partonic level. The strangeness suppression factor could then 
be a consequence of the large strange quark mass, compared to that of the light u and d 
quarks. In a partonic medium of temperature T > Th, the ratio of Boltzmann weights, 

7s ~ exp{-(ms - m„)/T} (31) 

would indeed lead to values around 0.5 for partonic temperatures around 200 MeV. Eq. 
(|3ll) implies that with increasing T, i.e., with increasing energy density of the excited 
vacuum, 7<. should approach unity. One way to check this is to see if the relative abundance 
of strange particles increases when one increases either the collision energy or considers 
collision fluctuations of higher than average secondary multiplicity. In both cases, the 
K/ vr ratio indeed grows considerably |]42| , ^ , as seen in Fig. |ll|. 



To complete the discussion of nucleon-nucleon collisions, we specify a{s) in Eq. 
From p — p and p — p experiments over a large range of collision energies, up to ^/s = 1 



TeV, one finds in the central rapidity region for the sum over all secondaries (which 
means predominantly pions). 



'dN, 



\ dy 

so that a{s) ~ \n{y/s/2m); here m denotes the nucleon mass. 



~ ln(v^/2m), (32) 

y~0 
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Figure 11: The K/ vr ratio (a) mp — p and p — p collisions as function of the cms collision 
energy and (b) 'm. p — p collisions at y/s = 1.8 TeV as function of the associated 
multiplicity [||]. 



After this brief survey of nucleon-nucleon interactions, we now return to nucleus- 
nucleus collisions in more detail. At high energies, each of the two incoming nuclei (for 
simplicity, we assume them to have equal mass numbers A) will appear Lorentz-contracted 
to a longitudinal thickness of 

dA = 2Ra i-^j (33) 



in the overall center of mass. Here s = 2m? + 2myP^^+rri? denotes the squared nucleon- 
nucleon collision energy, P the incident lab momentum per nucleon, and Ra — 1.12 A^^^ 
fm the nuclear radius. For c.m.s. collision energies of 20 GeV or more, (Ia is thus one 
fermi or less, so that the different nucleon-nucleon collisions in the course of the nucleus- 
nucleus collision essentially occur superimposed in the same space-time region. Consider a 
central A — A collision; in the process, each of the A projectile nucleons will interact one or 
more times with nucleons from the target nucleus. When the nuclei have passed through 
each other, we are thus left with 2 A 'wounded' nucleons, each of which eventually emits 



secondary hadrons |^5|. At present energies it is found that the number of secondaries 
is proportional to the number of wounded nucleons; it is not affected by the number of 
collisions which the nucleon has experienced Q. Since a proton-proton collision results in 
two wounded nucleons, we obtain from Eq. (|32D that the number of secondary hadrons in 
a central high energy A — A collision is given by 

(dNh\ ^ Ahi{^s/2m). (34) 



dy 



y~0 



^Note that this holds for 'soft' secondaries; the emission of 'hard' products, such as Drell-Yan dileptons 
or charmonium states, does depend on the number of coUisions. 
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In effect, this result is modified already in p— A collisions, since the nuclear target stops 
the projectile nucleon more than a nucleonic target does. This is taken into account by 
the form 



dy 



~ A"ln(v^/2m), (35) 

y~0 



with a ~ 1.15 for secondaries close to target rapidity in p— A data The multiplicity 



predicted by Eq. (|35D agrees quite well with results from S — S to Pb — Pb collisions at 
the CERN-SPS ||, H. 



3.2 Collision Evolution and Initial Conditions 

Statistical QCD, as described in Section 2, provides predictions for time-independent 
systems; we want to apply these to the rapidly evolving systems produced in nuclear 
collisions. Hence we have to consider in some detail the different evolution stages of such 
collisions. In the period of nuclear passage, which at high energies lasts a fermi or less, 
there is multiple scattering among the quark and gluon constituents of the colliding nu- 
cleons, leading to a rapid generation of entropy and thereby eventually to thermalisation. 
The separating nucleons thus leave behind a more or less equilibrated expanding medium 
which in its 'hot' early stages we expect to be a quark-gluon plasma. It cools off, under- 
goes the confining quark-hadron transition and after that forms hadronic matter, which 
finally freezes out into the observed hadronic secondaries. The aim of this section is to 
estimate the energy density of the medium in its early stage. 

Thermalisation is presumably attained by multiple scattering of the incident con- 
stituents and the subsequently produced secondaries. At high energies and early times, 
these constituents are quarks and gluons, so that the tool to study the resulting entropy 
production must be some form of parton cascade model based on perturbative QCD 
pO| , pl[| . The partons contained in the incident nucleons collide, produce further partons, 
and this cascade continues and eventually leads to a thermal medium. Such partonic 
cascade models allow the calculation of entropy production as a function of the local time 
after the collision, and the systems produced in high energy collisions appear to reach 
thermalisation after about one fermi. We shall return to parton cascade models in section 
6.3. 

To estimate the energy density, one has to model the expansion process. The sim- 



plest model 1 35] assumes that the finally observed hadrons have emerged from an initial 
interaction volume in free flow, i.e., without any collective motion. The resulting initial 
energy density in an A — A collision is 

,„ = ^f^V", (36) 



where po denotes the average (local) energy of the emitted secondaries and (dNh/dy) 
their multiplicity as given by Eq. (^3|). The initial interaction volume is determined by 
the transverse size Ra and the longitudinal extension as obtained from the thermalisation 
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time To, which is generally taken to be about one fermi. Combining this with Eq. 
we get 

eo ~ ln(yi/2m) [GeV/W] (37) 

for the initial energy density in a central A — A collision. Since pions are the dominant 
secondaries and have an average transverse momentum of about 0.35 GeV over a wide 
range of incident energies, we have po ~ 0.5 GeV. The estimates in Eq's. (^) and (|37D 
give averages over the transverse nuclear area, without taking into account any details of 
the nuclear distribution. In the Glauber formalism, the correct nuclear geometry can be 
taken fully into account (see e.g. |S^). For central A — A collisions, this results in little 
change, but for non-central or A — B collisions with A ^ B, the actual geometry becomes 
important. 

Before looking at the attainable energy densities, let us note where nucleus-nucleus 
collision experiments reaching the range of interest for our purposes can actually be per- 
formed. The only laboratories providing sufficiently energetic nuclear beams are the 
Brookhaven National Laboratory (BNL) near New York (USA) and the European Orga- 
nization for Nuclear Research (CERN) in Geneva (Switzerland). Both began experimenta- 
tion in 1986, using existing accelerators. BNL had the Alternating Gradient Synchrotron 
(AGS), designed for 30 GeV/c proton beams, CERN the Super Proton Synchrotron (SPS) 
for 450 GeV/c protons. The injectors available at that time allowed only the acceleration 
of nuclei containing equal numbers of protons and neutrons (A=2Z), so that the beams 
were restricted to light ions {A<40). Both laboratories have in the meantime built new 
injectors, allowing the acceleration of arbitrarily heavy nuclei in AGS and SPS. Moreover, 
both laboratories are presently constructing colliding beam accelerators which will bring 
them to very much higher collision energies. At BNL, the Relativistic Heavy Ion Collider 
(RHIC) is scheduled to start operating in 2000; at CERN, the Large Hadron Collider 
(LHC) is planned for some five years later. 

The available and planned facilities are summarized in Table 1. We have there also 
included the corresponding energy density values; they were obtained by a Glauber anal- 



ysis taking into account the correct nuclear geometry Note that the listed values 

are averages over the nuclear profile; the center of the interaction region can reach energy 
densities up to 30 % higher. The results using light ion beams are given for heavy targets 
[A ~ 200), the others for A — A with A ~ 200. The ln^/s dependence of the energy 
density (Eq. makes it rather difficult to achieve significant changes by varying the 
beam energy. The use of lighter nuclei is clearly more efficient. In Fig. |12|, we show the 
functional behaviour of e for A = 30 and A = 200, to give an indication of the ranges 
achievable at the different facilities, as well as their relation to the expected deconfinement 
threshold. 



3.3 The Onset of Deconfinement 



At the beginning of a high energy nucleus-nucleus collision, there are hard interac- 
tions between the partonic constituents of the participating nucleons. These partons are 
confined: their momentum distribution is that of partons within a nucleon, as determined 
e.g. in deep inelastic scattering. In particular, this means that if one quark in the nucleon 
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Machine 


Start 


lype 


Beam 


[GeV/A] 


AB 
^0 

[GeV/fm^] 


BNL - AGS 
CERN - SPS 


1986 
1986 


Fixed Target 
Fixed Target 




5 

19 


0.7 
1.6 


BNL - AGS 
CERN - SPS 


1992 
1994 


Fixed Target 
Fixed Target 


208p^ 


5 

17 


1.5 
3.7 


BNL - RHIC 
CERN - LHC 


2000 
~ 2005 


Collider 
Collider 


208p^ 


200 
5000 


7.6 
13 



Table 1: Experimental facilities for high energy nuclear collisions; the hght ion beam 
results are for heavy {A — 200) targets, the others for symmetric {A — A) collisions. 
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Figure 12: The variation of the initial energy density with incident collision energy for 
different experimental facilities. 



scatters, emitting (or absorbing) a gluon, others must follow the scattered quark to re- 
tain colour neutrality (Fig. |l^). As a result, high momenta k for the emitted gluon are 
suppressed, with a momentum distribution of the form 

x^(x) ~ (l-x)^+", (38) 

where x = kj ^ is the fractional momentum of the gluon and n denotes the number of 
other quarks or antiquarks which must follow the scattered one. 




Figure 13: Scattering of quarks confined to a nucleon. 

Deconfinement starts when the given environment leads to partons which are no longer 
subject to the confinement constraints just mentioned; with n = 0, this results in a 
hardening of the parton momentum distribution. The first step towards deconfinement is 
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thus the overlap of two hard nucleon-nucleon colhsions. If there is a secondary interaction 
between the two interacting pairs before the nucleons in the primary interactions have 
had the time to 'regroup' into physical nucleons, a hard gluon exchanged between the two 
pairs is not subject to confinement. An 'exogamous' system of this kind is therefore the 
first step towards a deconfined medium. 




Figure 14: 'Exogamous' interaction between two colliding nucleon systems (Pi,Ti) and 

(P2,T2). 

The interaction of two nucleons leads to the emission or absorption of gluons of trans- 
verse momenta fc^, with a corresponding transverse area of radius rj' ~ l/kx- In in- 
teractions of the type shown in Fig. |T^, partons from different collisions overlap in the 
transverse plane and can thus no longer be associated to a given nucleon. To form a con- 
nected medium of such deconfined partons, their density in the transverse plane has to be 
sufficiently high. This density increases when the mass number A (more collisions) and/or 
the collision energy (more partons) is increased; hence we need high energy collisions of 
heavy nuclei for the formation of a quark-gluon plasma. 



4. PROBES OF PRIMORDIAL MATTER 
4.1 Evolution 

We now assume that for sufficiently high energies, the deconfined partonic system 
quickly attains equilibrium; this is in accord with the parton cascade models mentioned 



above |5y, Subsequently, it cools off and then hadronizes. For an infinite medium, the 
actual transition process depends considerably on the order of the confinement transition, 
which so far remains uncertain. The systems produced in nuclear collisions are rather 
small, however, and in addition it is unlikely that any measurement will fix thermal 
observables precisely. Hence the effect of the transition will always appear in somewhat 
smeared-out form. On the hadronic side, it is possible for the medium to exist for some 
time as a system of interacting hadrons, which eventually freezes out. However, it is 
known that an ideal resonance gas leads to critical behaviour because of the composition 
law for resonances PBI, B^, and that an ideal gas of extended hadrons does so for geometric 



reasons 



53[| . Hence it is quite conceivable that the confinement transition leads directly 
into an ideal resonance gas. 

It is evident that the different evolution stages, from primary collisions to freeze-out, 
require very different probes for their investigation. Consider the two extremes. Pions 
appear at the very end, either directly or as decay products of resonances; hence they 
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can hardly tell us much about the very early phases of the system, which existed and 
disappeared long before pions and other large 'soft' hadrons were formed. Hard Drell-Yan 
dilcptons arc produced essentially at the time of the collision and then remain unaffected 
by the subsequent history of the system; hence they are not able to provide information 
on thermalisation and freeze-out. 

The basis for a quark-gluon plasma is the high parton density in the early stages after 
the collision; this effectively screens all confining long range forces, so that only short 
range interactions remain. Any deconfinement probe must therefore 

• be hard enough to resolve sub-hadronic scales, 

• be able to distinguish confinement and deconfinement, 

• be present in the early stages of the evolution, and 

• retain the information throughout the subsequent evolution. 

The last point requires that the probe should not be in thermal equilibrium with the later 
evolution stages, since this would lead to a loss of memory of the previous stages. 

So far, two types of probes satisfying these conditions fully or in part have been 
considered. 

• External probes are produced essentially by primary parton collisions, before the 
existence of any medium, and then indicate by their observed behaviour whether the 
subsequent medium was deconfined or not. The production of quarkonium states 
{J/ip, T) provides the best known example of such a probe. Short distance QCD 
predicts that the dissociation of quarkonia is possible only in a deconfined medium, 
since it requires harder gluons than those present in hadrons. Another example, 
quite hkely to gain considerable importance when the high energy coUiders RHIC 
and LHC come into operation, is the energy loss or attenuation of hard jets, which 
is expected to increase considerably in a deconfined medium. 

• Internal probes are produced by the quark-gluon plasma itself. Since they must leave 
the medium without being affected by its subsequent states, they should undergo 
only weak or electromagnetic interactions after their formation; the main candidates 
here are thermal dileptons or photons. However, since both can also be produced 

by a confined medium, they can serve as thermometer of the medium rather than as 
probe for its confinement status. So far, another problem is their identification: both 
dileptons and photons are produced abundantly in hadron decays, and separating 
out a possible thermal component appears very difficult. 

At the present stage of the investigation, quarkonium dissociation and jet quenching 
therefore appear the most promising direct signatures for deconfinement. 

Before concentrating on these hard deconfinement probes, we will first consider elec- 
tromagnetic probes in somewhat more detail. Real or virtual photons emitted during the 
evolution of the collision subsequently undergo no (strong) interactions with the medium 
and hence refiect its state at the time they were produced. On the other hand, they are 
emitted during the entire collision evolution, and by different dynamical mechanisms at 
different stages: 
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• Early hard parton interactions produce hard photons and Drell-Yan dileptons; these 
provide information about the initial (primary or pre-equilibrium) stages. 

• Thermal photon and dilepton emission by the medium, through quark or hadron 
interactions, occur through its entire evolution, and hence give information about 
the successive stages, from quark-gluon plasma to final hadronic freeze-out. 

• Hadrons produced at any point of the hadronic stage, from the quark-hadron tran- 
sition to freeze-out, can decay and thereby emit photons or dileptons; depending on 
the hadron decay time, they provide information about dense interacting hadronic 
matter or about hadrosynthesis at the end of the strong interaction era. 

Each of these formation mechanisms provides information about specific aspects of the 
different collision stages. 

Drell-Yan dileptons and hard photons are the tools to study the effective initial state 
parton distributions; in particular, they will show any nuclear modifications (shadowing, 
anti-shadowing, coherence effects) of these distributions. They also indicate the initial 
state energy loss and the initial state px broadening suffered by partons in normal nuclear 
matter. Since they do not undergo any final state strong interactions, they moreover 
provide a reference for the effect of the produced medium on quarkonium states or jets. 

Thermal emission can in principle be used to determine the temperature at the different 
evolution stages. The functional form of thermal spectra, 

dN/dk^ ~ e-''-'/^ (39) 

for photon momenta, or the corresponding distributions in the dilepton mass Mi+i-, indi- 
cate the temperature T of the medium at the time the signal was emitted. Since the form 
(|39|) for thermal production is the same for a hadronic medium and for a quark-gluon 
plasma, it cannot specify the nature of the emitting medium, only its temperature; this 
feature finds support in various cascade model calculations 0, The crucial thermo- 
metric problem is to find a "thermal window" , since the measured spectra are dominated 
at high photon momenta or dilepton masses by hard primary reactions and at low mo- 
menta or masses by hadron decay products. As mentioned, thermal photon or dilepton 
emission has so far not been identified, perhaps because of the dominance of the hadronic 
stage at present energies. 

The dileptons produced by hadron decays can probe hadronic in-medium modifica- 
tions, provided the hadrons actually decay within the medium. The p, with a half-life of 
about a fermi, appears to be the best candidate for such studies. Chiral symmetry restora- 
tion is expected to change the properties of hadrons as the temperature of the medium 
approaches the restoration point |Q; hence such in-medium changes are of particular 
interest, since they might be the only experimental tool to address the chiral aspects of 
deconfinement. If at the onset of chiral symmetry restoration, the mass of the p decreases 



sufficiently much [55 



""^^^^ ' as T ^ r„ (40) 



m,(r = 0) 

then this should lead to an observable enhancement of the low mass dilepton spectrum. 



There are first indications for such an enhancement |Q - |]58[ ; however, so far alternative 
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explanations based on in-medium hadron interactions, leading to a change of resonance 
widths, can also account for the effect [^. Moreover, a temperature dependence of the 
form ( ^OD has so far not been corroborated by lattice studies of hadron masses |Q . 

In the following two subsections, we shall discuss how the behaviour of quarkonium 
states and of hard jets can provide information about the confinement status of a given 
medium. In order to apply the results as probes in actual nuclear collision experiments, 
the production of the probes in the collision process must be understood as well, which 
clearly complicates matters. In section 5, we shall carry out a full analysis of this type for 
charmonium production and suppression in nuclear collisions; for this, there now exists a 
wealth of beautiful data IBlll. 



4.2 Quarkonium Dissociation 

We begin our study of deconfinement probes by restating the basic question: Given a 
box of unidentified matter, determine if the quarks and gluons which make up this matter 
are in a confined or deconfined state. It is thus not evidence for quarks and gluons that we 
look for, but rather the confinement/deconfinement status of these elementary building 
blocks of all forms of matter. 

As quarkonium prototype, consider the J /ip; it is the 15* bound state of a charm quark 
(m-c ~ 1.4 GeV) and its antiquark, with Mj/^ ^3.1 GeV. Its usefulness as deconfinement 
probe is easily seen. If a J/'?/' is placed into a hot medium of deconfined quarks and gluons, 
colour screening will dissolve the binding, so that the c and the c separate. When the 
medium cools down to the confinement transition point, they will therefore in general be 
too far apart to see each other. Since thermal production of further cc pairs is negligibly 
small because of the high charm quark mass, the c must combine with a light antiquark 
to form a and the c with a light quark for a D. The presence of a quark- gluon plasma 
will thus lead to a suppression of J / ip production . 

We shall now first consider this J/tp suppression in terms of colour screening, i.e., as 
consequence of global features of the medium, and then turn to a microscopic approach, 
in which the bound state is assumed to be broken up by collisions with constituents of 
the medium. 

4.2a. J/ip Suppression by Colour Screening 

Because of the large charm quark mass, the charmonium spectrum can be calculated 
with good precision by means of the Schrodinger equation 

[2m, + — + Vir)]^n,i = M„,^„,z, (41) 

where the potential V{r) = or — ajr contains a confining long-distance part or and 
a Coulomb-like short- distance term ajr. For different values of the principal quantum 
number n and the orbital quantum number /, the masses M^^i and the wave functions 
^n,i(^) of different charmonium states J/"?/", ip' ., ... in vacuum are given in terms of the 
constants mc, a and a. 
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In a medium, the potential becomes screened, as we had seen above in Eq. we thus 
obtain 

(J ol 



\/(r,/i) = -[l 



a 

—e 
r 



(42) 



where is the screening mass and r^j = the 'Debye' colour screening radius. Screening 
is a global feature of the medium, shortening the range of the binding potential. Once 
11 becomes sufficiently large, the bound states begin to disappear, starting with the most 
weakly bound; hence for // > /x^, the bound state i is no longer possible [ p^ . 
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^'(2S) 


Xc(lP) 


T(1S) 


T'(2S) 


T"(3S) 


X6(1P) 


X', (2P) 


^id [GeV] 


0.68 


0.35 


0.35 


1.52 


0.64 


0.43 


0.59 


0.42 




1.2 


1 


1 


2.7 


1.1 


1 


1.1 


1 


AE [GeV] 


0.64 


0.06 


0.24 


1.10 


0.54 


0.20 


0.67 


0.31 



Table 2: Quarkonium Dissociation by Colour Screening 

With the help of finite temperature lattice QCD, we can now determine the screening 
mass niT) as function of the temperature T. With the approximate relation fi ^ 3.3 T 
and Tc ~ 0.17 GeV, one obtains the melting pattern for the most important cc and bb 
states summarized in Table 2; we see that while both ip' and Xc melt at the critical 
deconfinement point, the J/ip, being smaller and more tightly bound, survives to about 
2Tc and hence about twice the critical energy density. A similar pattern holds for the 
bottonium states; the dissociation values of the screening masses for Xb T" have not 
yet been determined in potential theory. With increasing temperature, a hot medium will 
thus lead to successive quarkonium melting, so that the suppression or survival of specific 
quarkonium states serves as a thermometer for the medium, in much the same way as the 
relative intensity of spectral lines in stellar interiors measure the temperature of stellar 
matter ||6^ . Note, however, that other possible sources of quarkonium dissociation have to 
be considered before the method becomes unambiguous. This problem will be addressed 
in detail in section 5. 

4.2b. J/^ Suppression by Gluon Dissociation 

The binding energy of the J/ip, i.e., the energy difference between the J/ip mass and 
the open charm threshold, AEj/^ = 2Md — Mjj^ ~ 0.64 GeV is considerably larger than 
the typical non-perturbative hadronic scale Aqcd — 0.2 GeV. As a consequence, the size 
of the J /tp is much smaller than that of a typical hadron, rjj^ ~ 0.2 fm -C Aqqj-, = 1 fm. 
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Figure 15: J/ip dissociation by hadron interaction. 



Hence the J /if) is a. hadron with characteristic short-distance features; in particular, rather 
hard gluons are necessary to resolve or dissociate it, making such a dissociation accessible 
to perturbative calculations. J/ip collisions with ordinary hadrons made up of the usual 
M, d and s quarks thus probe the local partonic structure of these 'light' hadrons, not their 
global hadronic aspects, such as mass or size. It is for this reason that J/ip^s can be used 
as a confinement /deconfinement probe. 

This can be illustrated by a simple example. Consider an ideal pion gas as a confined 
medium. The momentum spectrum of pions has the Boltzmann form f{p) ~ exp — {\p\/T), 
giving the pions an average momentum = 3 T. With the pionic gluon distribution 
function xg{x) ~ (1 — x)^, where x = k/p denotes the fraction of the pion momentum 
carried by a gluon, the average momenta of gluons confined to pions becomes 

{\k\)conf ^ 0.6 T. (43) 

On the other hand, an ideal QGP as prototype of a deconfined medium gives the gluons 
themselves the Boltzmann distribution f{k) ~ exp — {\k\/T) and hence average momenta 

(|A:|)deconf = 3T. (44) 

Deconfinement thus results in a hardening of the gluon momentum distribution. More 
generally speaking, the onset of deconfinement will lead to parton distribution functions 
which are different from those in vacuum, as determined by deep inelastic scattering 
experiments. Since hard gluons are needed to resolve and dissociate J/ip's, one can use 
J/^/^'s to probe the in-medium gluon hardness and hence the confinement status of the 
medium. 

These qualitative considerations can be put on a solid theoretical basis provided by 
short- distance QCD |^ - [^. In Fig. |l^ we show the relevant diagram for the calculation 



of the inelastic J/ip-hadion cross section, as obtained in the operator product expansion 
(essentially a multipole expansion for the charmonium quark- ant iquark system). The 
upper part of the figure shows J/ ip dissociation by gluon interaction; the cross section for 
this process, 

(Tg.j/^ ~ (A; - /\E^f'^k-\ (45) 

constitutes the QCD analogue of the photo-effect. Convoluting the J /tp gluon-dissociation 
with the gluon distribution in the incident hadron, xg{x) ~ 0.5(1 — x)^^"^, we obtain 

cTh-j/^ ^ ageom(l - Ao/A)"+=^-^ (46) 



28 



1e-01 













/'n - \|/ 









k [GeV] 

Figure 16: J/ip dissociation by gluons and by pions; k denotes the momentum of the 
projectile incident on a stationary J/ip. 




Figure 17: Elastic J/ip photoproduction data, compared to the prediction obtained from 
short-distance QCD and VDM |TI1. 
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for the inelastic J/ip-hadron cross section, with A ~ {s — MV)/M.^ and Aq — (M/i + Ai?^); s 
denotes the squared J/'i/'-hadron coUision energy. In Eq. (^q) , Ugeom — const, ~ 2—3 mb 
is the geometric cross section attained at high colhsion energies with the mentioned gluon 
distribution. In the threshold region and for relatively low collision energies, ah-j/ij) is very 
strongly damped because of the suppression (1 — xY^^ of hard gluons in hadrons, which 
leads to the factor (1 — Aq/A)"^^'^ in Eq. (^). In Fig. |16|, we compare the cross sections for 
J/ip dissociation by gluons ( "gluo-effect" ) and by pions = 2), as given by Eq's ( ^5D and 
(|46|). Gluon dissociation shows the typical photo-effect form, vanishing until the gluon 
momentum k passes the binding energy AE^; it peaks just a little later and then vanishes 
again when sufficiently hard gluons just pass through the much larger charmonium bound 
states. In contrast, the J/ip-hadron cross section remains negligibly small until rather 
high hadron momenta (3-4 GeV). In a thermal medium, such momenta correspond to 
temperatures of more than one GeV. Hence confined media in the temperature range of 
a few hundred MeV are essentially transparent to J/ip^s, while deconfined media of the 
same temperatures very effectively dissociate them and thus are J/V'-opaque. 

The formalism leading to Eq. (HB) is also applicable to J/ip photoproduction, where 



it can be compared directly to data |]70[. An alternative relation between J/ip photopro- 
duction and J/ip break-up by hadrons is obtained through a combination of the vector 
meson dominance model and the short- distance QCD approach [^]. In both cases one 



finds clear evidence for the strong threshold suppression due to the form of the gluon 
distribution function; an illustration is given in Fig. |T^. In contrast, several phenomeno- 
logical models based on effective hadronic Lagrangians |]72| - suggest a cTg^j/^ which 
is very much larger than Eq. (|l6|) in the threshold region. It seems unlikely that such 
models can survive a comparison to J/ip photoproduction. 

The situation for Xc's is quite similar, except for an earlier gluon dissociation threshold 
due to the lower Xc binding energy of about 250 MeV. Corresponding differences occur 
for the various bottonium states. The difference in binding energies thus provides the 
microscopic basis for the successive melting of different quarkonium states noted above. 
The binding energies for the quarkonium states considered above are included in Table 2. 
We note in particular that for the AE is almost negligible (about 40 MeV), so that 
there cannot really be any difference in its dissociation by confined or deconfined media. 
In other words, any form of strongly interacting matter is expected to be V^'-opaque. 

We can thus define a schematic test to determine the confinement status of a box 
of unidentified matter. We first shine a ip' beam at it: if it is transparent to this, the 
box does not contain strongly interacting matter; otherwise it does. In that case we 
repeat the procedure with a J/ip beam: if this is unaffected, the medium in the box is 
confined; if the beam is attenuated, it is deconfined. The problem of an unambiguous 
deconfinement test is thus in principle solved: there is J/ijj dissociation if and only if the 
medium is deconfined. However, the test pre-supposes the existence of a prepared strongly 
interacting medium and the availability of J/ijj and ip' beams as probes, and in nuclear 
interactions, both have to be prepared by the collision process. The actual application of 
such a test is therefore more complex; it will be taken up in detail in section 5. 
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(a) (b) 

Figure 18: Some lowest order diagrams for jet production in hadron collisions. 
4.3 Jet Quenching 

One important aspect in the use of quarkonium states as deconfinement test is that in 
a collision with usual 'big light' hadrons, the small quarkonia probe only the local partonic 
structure of their collision partner, not global features such as mass or size. This aspect 
is also the basis for using hard jets to study the confinement status of a given medium. 

Hard jets are produced through hard partonic interactions in the very early stages of 
nucleon-nucleon collisions; some typical processes are illustrated in Fig. |18|. In nuclear 
collisions, the produced gluons or quarks leave their formation point with a very high 
momentum transverse to the collision axis, and in doing so, pass through the primary 
nuclear matter and through whatever secondary medium is produced by the collision. 
The effect of such ge through a quark-gluon plasma has in recent years been 



studied in considerable detail, leading to rather well-defined predictions - [[f^ . 

An electric charge, passing through matter containing other bound or unbound charges, 
loses energy by scattering. For charges of low incident energy E, the energy loss is largely 
due to ionization of the target matter. For sufficiently high energies, the incident charge 
scatters directly on the charges in matter and as a result radiates photons of average 
energy uj E. Per unit length of matter, the 'radiative' energy loss due to successive 
scatterings, 

dE ^ 

E 47 

dz 

is thus proportional to the incident energy. 

This probabilistic picture of independent successive scatterings breaks down at very 
high incident energies |8^. The squared amplitude for n scatterings now no longer fac- 
torizes into n interactions; instead, there is destructive interference, which for a regular 
medium (crystal) leads to a complete cancellation of all photon emission except for the 
first and last of the n photons. This Landau-Pomeranchuk-Migdal (LPM) effect greatly 
reduces the radiative energy loss. 

The physics of the LPM effect is clearly relevant in calculating the energy loss for 
fast colour charges in QCD media. These media are not regular crystals, so that the 
cancellation becomes only partial. Let us consider the effect here in a heuristic fashion; 
for details of the actual calculations, see |7^, |7^. The time tc needed for the emission of 



a gluon after the scattering of a quark (see Fig. |T9D is given by 



I E E 
Jp2jpi 2P'k ^ ' 
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Figure 19: Gluon emission after scattering. 
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Figure 20: Energy loss in incoherent and coherent interactions. 
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Figure 21: Energy loss in coherent interactions as function of the thickness L of the 
medium 
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in the rest system of the scattering center, where measures how far the intermediate 
quark state is off-shell; on-shell quarks and gluons are assumed to be massless, and E/^/P^ 
is the 7-factor between the lab frame and the proper frame of the intermediate quark. 
For gluons with ki » kr, we thus get 

tc ^ ^- (49) 

If the passing colour charge can interact with several scattering centers during the for- 
mation time of a gluon, the corresponding amplitudes interfere destructively, so that in 
effect after the passage of n centers over the coherence length Zc-, only one gluon is emitted, 
in contrast to the emission of n gluons in the incoherent regime. Nevertheless, in both 
cases each scattering leads to a fc^-kick of the charge, so that after a random walk past n 
centers, k"^ ~ n. Hence 

4 ^ /i'y, (50) 

where A is the mean free path of the charge in the medium, so that A > 1 counts the 
number of scatterings. At each scattering, the transverse kick received is measured by the 
mass of the gluon exchanged between the charge and the scattering center, i.e., by the 
screening mass /i of the medium. From Eq. ^ we have 

UJ 

Zc ^ p , (51) 
so that the formation length in a medium characterized by /z and A becomes 

(52) 

For the validity of Eq. (0), the mean free path has to be larger than the interaction range 
of the centers, i.e., A > 

The energy loss of the passing colour charge is now determined by the relative scales of 
the process. If A > Zc, we have incoherence, while for X < Zc there is coherent scattering 
with destructive interference. In both cases, we have assumed that the thickness L of the 
medium is larger than all other scales. When the coherence length reaches the size of the 
system, z^ = L, effectively only one gluon can be emitted. This defines a critical thickness 
L^{E) = (EA//i2)V2 at fixed incident energy E, or equivalently a critical E^ = ^Jt^L"^ / A for 
fixed thickness L; for L > L^, there is bulk LPM-behaviour, below there are finite-size 
corrections. 

We are thus left with three regimes for radiative energy loss. In case of incoherence, 
Zc < there is the classical radiative loss 

(53) 

dz TT A 

where is the strong coupling. In the coherent region, A > z^-, the energy loss is given 
by the LPM bulk expression when L > L^. 

Q^, I ,.2 TP 

(54) 
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The resulting reduction in the radiative energy loss dE/dz is illustrated in Fig. 20. Note 
that in earlier estimates the energy loss due to interactions of the gluon cloud accompa- 
nying the passing colour charge had been neglected |]7^; this led to a considerably smaller 
energy loss, proportional to In E instead of \fE. Finally, in a medium of thickness L < Lc, 
there is less scattering and hence still less energy loss. Eq. (E^ can be rewritten as 



--T- —^L,iE), (55) 
az IT A 

dE 3asfi'\ 
- — ~ -L 56 

dZ TT A 

as the energy loss in finite size media with L < L^.. The resulting variation of the radiative 



and for L < L^, this leads to 

,2 



energy loss with the thickness of the medium is shown in Fig. 21, with saturated (i.e., 
bulk) LPM behaviour setting in for L > L^. 

Eq. (p6D has been used to compare the energy loss in a deconfined medium of tem- 



perature T = 0.25 GeV to that in cold nuclear matter of standard density |81]. For the 



traversal of a medium of 10 fm thickness, estimates give for the total energy loss 

rlQ fm 

AE= / dz^ (57) 

JO fm dz 

in a quark-gluon plasma 

-AE,gp~30 GeV, (58) 
corresponding to an average loss of 3 GeV/fm. In contrast, cold nuclear matter leads to 

-AE,„„^2 GeV (59) 

and hence an average loss of 0.2 GeV/fm. A deconfined medium thus leads to a very much 
higher rate of jet quenching than confined hadronic matter, as had in fact been suggested 
quite some time ago p2| . 

To convert this into an operational signature for deconfinement, considerably more 



work is needed [p3| , p4| . Some further basic ingredients needed for this effort will become 



clear in the next section, where the corresponding task is carried out for J/-ip suppression. 



5. J/^ SUPPRESSION IN NUCLEAR COLLISIONS 

In this section, we shall present a realistic deconfinement study, based on the analysis of 
J/ ijj production in nuclear collisions. On one hand, this provides a nice application of the 
general ideas presented above; on the other hand, it illustrates the steps necessary in any 
other, future deconfinement study. First, the production of the probe must be understood 
in hadronic collisions, both theoretically and experimentally. We therefore begin with the 
hadropro duct ion of charmonium. Next, one has to understand the modifications arising 
when the production occurs in a confined medium, for which p — A collisions provide the 
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experimental reference. With the probe thus (a) prepared and (b) gauged to confined 
matter, it can be apphed in an environment in which there might be deconfinement. 



5.1 The Hadroproduction of Charmonium 

The first stage of charmonium formation in hadronic colhsions is the production of a 
cc pair; because of the large quark mass, this process can be described by perturbative 
QCD (Fig. |2^). A parton from the projectile interacts with one from the target; the 
(non-perturbative) parton distributions within the hadrons are determined empirically in 
other reactions, e.g. by deep inelastic lepton-hadron scattering. Initially, the cc pair will 
in general be in a colour octet state. It subsequently neutrahses its colour and binds to a 
physical resonance, such as J/V', Xc or '4'' ■ Colour neutralisation occurs by interaction with 
the surrounding colour field; this and the subsequent resonance binding are presumably 
of non-perturbative nature. The colour evaporation model ||8^ provides the simplest and 



most general phenomenological approach to colour neutralisation. In the evaporation 
process, the cc pair can either combine with light quarks to form open charm mesons {D 
and D) or bind with each other to form a charmonium state. 

The basic quantity in this description is the total sub-threshold charm cross section, 
obtained by integrating the perturbative cc production cross section over the mass interval 
from 2mc to 2m£). At high energy, the dominant part of acc comes from gluon fusion (Fig. 



2^a), so that we have 

2m,£) 



2mc 



ds / dxidx2 Qpixi) gt{x2) cr{s) 5{s — xiX2s), (60) 



with gp{x) and gt{x) denoting the gluon densities, Xp and Xt the fractional momenta 
of gluons from projectile and target, respectively; a is the gg cc cross section. In 
pion-nucleon collisions, there are also significant quark-antiquark contributions (Fig. ^c), 
which become dominant at low energies. The essential claim of the colour evaporation 
model is that the production cross section of any charmonium state i is given by 

c^i(s) = fi (Tccis), (61) 

where fi is an energy-independent constant to be determined empirically. It follows that 
the energy dependence of any charmonium production cross section is predicted to be 
that of the perturbatively calculated sub-threshold charm cross section. As a further 
consequence, the production ratios of different charmonium states 

^ = A = const. (62) 

must be energy- independent. We note that only a small part of the total subthreshold 
cross section acc goes into charmonium formation. In accord with perturbative open 
charm calculations, the remainder (more than 90 %) leads to DD production; the missing 
energy needed to reach the DD threshold is obtained by interaction with the colour field. 
The cross section for open charm production is thus directly predicted by the form of Eq. 



35 




(a) 




(b) 




(c) 

Figure 22: Lowest order diagrams for cc production in hadronic collisions, through (a,b) 
gluon fusion and (c) quark-antiquark annihilation. 
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Figure 23: Open charm hadroproduction at different cms energies, compared to PYTHIA 
calculations [061. 
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(60), with the integration essentially running to a/s, instead of ending at 2m£). In Fig. 23 



we see that the result agrees well with experiment 

The predictions of the colour evaporation model have been compared in a comprehen- 
sive survey to the available data, using parton distribution functions which include 
the new HERA results |8^. In Figs. ^ and |2^, we see that the energy-dependence 
is well described for both J/ijj and T production; for J/ip production, the normalisation 
coeffficient is /j/^=0.025. The T results are obtained for the sum of T, T' and T" decay- 
ing into dimuons, with 5/x = 1.6 x 10~^ for the normalisation; here the branching ratios 
cannot be directly removed. In the fixed-target/ISR energy range, the results from the 
two different sets of parton distributions coincide; for the J/ip at LHC energies, there is 
some spread due to scale uncertainties in the parton distributions, which hopefully can be 
removed by more precise DIS data. For the T production, we have already now data up 
to 1.8 TeV, and in Fig. ^ they are seen to agree very well with the prediction obtained 
using the "low energy" value Bfr = 1.6 x 10~^. 



In Figs. |2g and the predicted energy-independence of production ratios is found to 
hold as well, again up to Tevatron energy. Here it should be noted that the CDF data for 
the ratio ip' /{J/ip) are taken at large transverse momenta (5 < < 15 GeV), while the 
lower energy data are integrated over p^, with the low px region dominant. Hence colour 
evaporation appears to proceed in the same way at both small and large pt- 

With the relative fractions fi for the different quarkonium states determined, colour 
evaporation also predicts the longitudinal and transverse momentum distributions for the 
hadroproduction of these states. In both cases, the agreement is quite good fS^, @ • 

The colour evaporation model thus provides a viable phenomenological description of 
the hadroproduction of quarkonia, leading to quantitative predictions up to the highest 
energies under consideration. However, it does not give a theoretical basis for the frac- 
tions fi of the hidden charm cross section, and hence one cannot use it to determine the 
space-time evolution of quarkonium production. For the production in nuclear collisions, 
this is crucial, as we shall see shortly, and so we have to consider more detailed dynam- 
ical mechanisms for hadroproduction. The first such attempt, the colour singlet model. 



assumed that the colour neutralisation of the cc is perturbatively calculable [QJ]. The re- 
sulting well-defined predictions were found to disagree by orders of magnitude with high 
energy data from FNAL |]92|, so that non-perturbative features seem to be essential for 



quarkonium production 94|. This has led to extensive studies of possible mechanisms; 



see for recent reviews. 

We shall here concentrate on one resulting proposal, the colour octet model p3| , |94 



keeping in mind, however, that only rather general features of such a model will enter the 
subsequent applications to nuclear collisions. The basis of the colour octet model is that 
any quark-antiquark bound state is in reality a superposition of the basic qq state and 
fluctuations involving in adddition one or more gluons ("higher Fock space components"). 
We can thus decompose the J /ip state 

= aol(cc)i) + ai\{cc)8g) + a2\{cc)igg) + a'2\{cc)8gg) + ... (63) 
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Figure 26: The ratio of J/ip to ip' hadropro duct ion as function of the cms energy [p7 
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into a pure cc colour singlet component (^5"!), into a component consisting of cc colour 
octet (^Sq or ^Pj) plus a gluon, and so on. The higher Fock space coefficients correspond 
to an expansion in the relative velocity v of the charm quarks. As we shall see below, this 
corresponds in coordinate space to an expansion in terms of components of decreasing 
spatial extension. For the wave function of the J/ip, the higher components correspond 
to (generally small) relativistic corrections. In high energy J/ip production, however, the 
higher Fock space components of charmonium states play an important role. The short 
time available before confinement constraints appear, favours colour neutralisation of the 
(cc)8 by gluons already present. Since (cc) production at high energies occurs at small 
X, the density of such comoving gluons is in fact high, and so the higher Fock space 
components become dominant. 

Analogous Fock space decompositions hold for the other charmonium states |^3|, ^ 



In all cases, the ffist higher state consists of a colour octet cc plus a gluon. For the ip', 
the next-to- leading terms again contain a colour octet (^Pj or ^5*0 cc) plus a gluon; for 
the x's, a ^5*1 colour octet cc is combined with a gluon. 

This sheds some light onto the unspecified colour evaporation process. When the 
colour octet cc leaves the field of the nucleon in which it was produced, it will in general 
neutralise its colour by combining non-perturbatively with an additional coUinear gluon, 
thus producing the {cc)sg component of the J/ip or the other charmonium states (Fig. 
P^). After a "relaxation time" rg, the {cc)sg will absorb the accompanying gluon to revert 
to the dominant (cc)i charmonium mode. Note that we are here considering those cc 
pairs which will later on form charmonia. The (cc)8 could as well neutralise its colour by 
combining with a light quark-antiquark pair, but this would result predominantly in open 
charm production. 
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Figure 28: J/ip production in p — p collisions. 

Charmonium production in hadronic collisions thus inherently involves different space- 
time scales [93|. The ffist step is the creation of a heavy cc pair, which takes place in a 
very short time, Tpcrt — l/(2mc). The colour octet cc state combines with an additional 
gluon to form a pre-resonance (cc)85' to neutralise its colour and yield a resonance state of 
the correct quantum numbers. The lifetime rg of this state is determined by the virtuality 
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of the intermediate cc state. In the rest frame of the cc, this is approximately p 



^ (64) 



where A = [{p + k)"^ — m^] = 2pk, with p denoting the four momentum of a quark in the 
final J/ip, k that of the accompanying gluon. Because of the large charm quark mass rric, 
the time uncertainty associated with gluon absorption becomes 

where ko is the gluon energy. In a confining medium, k > Aqcd, with Aqcd — 0.20 — 0.25 
GeV. We thus encounter a new scale (2mcAQCD)"^^^, with Aqcd ^ (SmcAqco)"'^^^ ^ 
it makes rg ^ Tport- For J/ip production at mid-rapidity of a nucleon-nucleon collision, 
the colour neutralisation time becomes 

t8~r8[l + (PA/2m,)2]V' (66) 

in the rest frame of target or projectile nucleon, with Pa denoting the momentum of the 
cc pair in this frame. As seen from the nucleon, colour neutralisation of fast cc pairs will 
thus take a long time. Equivalently, the (cc)85' travels in the time tg a long distance, 

ds ~ TsiPA/2m,), (67) 

in the rest frame of the nucleon. 

Moreover, the production of high energy charmonia implies the production of compos- 
ite and hence extensive states {ccg)i = |(cc)85')j; here the subscript i again refers to the 
specific charmonium state in question [J/ip, Xc '^')- Although the intrinsic transverse 
size of the {ccg)i depends in principle on the quantum numbers of the state i, confinement 
constraints on the gluon lead to a universal (ccg) size rg. Let us consider this in more 
detail. 

The formation process of the cc pair inside the (ccg) state (see Fig. |28|) makes this 
pair very compact, with a spatial extension of about l/2mc. The size of the (ccg) is thus 
essentially determined by the softness of the gluon, which is restricted hj k > Aqcd- The 
effect of this can be estimated by noting that in the time Tg obtained above, a gluon 
can propagate over a distance rg ~ (2mcAQCD)~^'^^- Equivalently, we can consider energy 
conservation in the passage from {ccg) to the leading Fock space component (cc)i; the 
non-relativistic form for heavy quarks leads to 

^ - h, (68) 



2mc 

where p is the quark three-momentum in the charmonium cms. /^From the lowest allowed 
gluon energy ko = Aqcd we thus obtain the intrinsic size 

r8 ~ - ~ ^=^= ~ 0.20 - 0.25 fm. (69) 
P J2mcAQCD 
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Since this size is determined only by the {ccj^g composition of the next-to- leading Fock 
space state, the compactness of the produced {cc)s in this state and the gluon momentum 
cut-off in confined systems, it is essentially the same for the different charmonium states. 
The Fock space state [ccg) thus constitutes something like a gluonic hard core present 
in all charmonium states. Its size is seen to be approximately that of the ground state 
charmonium J/ip, while the higher excited charmonium states are larger. The hard gluon 
core in these just corresponds to the fact that in extended bound states the emission (or 
absorption) of gluons with momenta bigger than Aqcd can occur only in a sufficiently 
small inner region. In other words: while the basic (cc)i state gives for different quantum 
numbers quite different spatial distributions, the much more localised higher Fock space 
states are of universal size (2mcAQCD)^^^^- In fhe large quark mass limit, however, the 
bound state radii become smaller than (2mQAQCD)~^^^, so that the universality would 
then be lost. 

In summary: charmonium hadroproduction is based on a hard perturbative process, 
such as gluon fusion into a cc pair. The subsequent colour neutralisation brings in non- 
perturbative features, leading to the formation of a pre-resonance state, such as the first 
higher Fock state component {ccjsQ- Confinement constraints result in a universal pre- 
resonance charmonium radius of some 0.25 fm and an associated life time of the same size. 
In this evolution era of charmonium production, the different resonant states {J/ip, Xc, 4'') 
thus appear essentially indistinguishable. This has striking consequences for charmonium 
production on nuclear targets. 

5.2 Pre- Resonance Suppression 

The colour evaporation model alone, without further dynamical input, predicts that 
the charmonium production rate in high energy p — A collisions is A times that in p — p 
collisions at the same energy, just as it is for open charm production. We assume here 
that the hard process is dominated by gluon fusion; hard processes involving quarks would 
lead to a small modification due to isopin effects. Experimentally, however, the presence 
of the nuclear target medium is known to reduce J/ ip production rates in p — A collisions 
significantly, up to about 40% relative to those m. p — p interactions P^'lHli -^^S- 



P9| we show the most recent data ||100|| . The suppression m. p — A collisions was the 
first indication that the presence of a strongly interacting medium influences the rate of 
J/ip production. The nuclear medium introduces space-time scales, and hence we have 
to consider the space-time structure involved in the production process. 

The mentioned experiments were carried out with incident protons with momenta 
between 200 and 800 GeV/c. This gives the nascent J/ip momenta of 30 GeV/c or 
more in the target rest frame. As a result, the transition ccg J/'^P, Xc or ip' occurs 
outside the target nucleus; the nuclear medium of the target sees only the passage of the 
corresponding pre-resonance states. Since these have essentially the same size and life- 
time for all charmonium states, the observed attenuation of the production rates should be 



the same for J/ijj as for ip'. This prediction is indeed found to be quite well satisfied ||105 
(see Fig. ^). Earlier attempts to explain charmonium suppression inp — A interactions in 
terms of the absorption of physical J/ ip states ||101|] had encountered difficulties precisely 



because of this feature. The equal attenuation of J/ip and ip' is a natural consequence 
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Figure 29: J /if) production in p — A vs. p — p collisions [|100|] . 



of pre-resonance absorption; it can never be obtained for the physical J/ip and ip' states 
with their very different geometric sizes. 

The cross section for the dissociation of the ccg state through collisions in nuclear 
matter can be estimated theoretically ||102|| ; it can also be determined directly from p — A 
data Since the components of the ccg are colour octets, the gluon interaction with 



this state will be 9/4 times stronger than the gluon interaction with the fundamental cc, 
whose components are colour triplets. Hence we expect the ccg dissociation cross section 
to be about 9/4 times the high energy J/ip dissociation cross section, i.e., around 6 mb. 

For an experimental determination, we note that a ccg pair formed at point zo in the 
target nucleus has a survival probability 



S. 



ccg 



exp 



dz pa{z) a, 



^0 



ccg 



-N 



(70) 



where the integration covers the path remaining from zq out of the nucleus. The traversed 
medium of nucleus A is parametrized through a Woods-Saxon density distribution pa{,z), 
and by comparing S^^ with data for different targets A, the dissociation cross section for 



ccg — N interactions is found to be ||103 



0", 



ccg- 



-N 



6.9 ± 1.1 mb. 



(71) 



In Fig. ^ it is seen that pre-resonance absorption with this cross section agrees well with 
all presently available p — A data on J / ip production. 

We now turn to nucleus-nucleus collisions; here there will certainly also be pre- 
resonance absorption in nuclear matter. However, in addition to the target and projectile 
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Figure 30: The relative 74-dependence of J/ip and ip' production [|105|| . 
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Figure 31: Schematic view of J/tp production in A — B colhsions. 

nuclei, there could now be a substantial amount of produced 'secondary' medium (Fig. 
^). Our ultimate aim is to look for such a medium and test its confinement status. 

The survival probability of a pre-resonance charmonium state in an A — B collision at 
impact parameter b is given by 

{roc roo | 

- I dz Pa{z) (Tccg-N g dz Pb{z) (Tccg-N \ , (72) 

in extension of Eq. (ffOl) . Here Zq specifies the formation point of the ccg within nucleus 
A, Zq its position in B. Since experiments cannot directly measure the impact parameter 
6, we have to specify how Eq. (ff^) can be applied to data. 

The Glauber formalism allows us to calculate the number N^^{b) of participant 
('wounded') nucleons for a given collision. The number of secondary hadrons produced in 
association with the observed J/ip is at present energies (up to SPS energy) found to be 



44 



1.4 



1.2 



coP 0.8 



C/D 



P 0.6 



0.4 
0.2 










L 


\ 

L 


\ 



"I WW 



{A p-A (450 GeV) 
O p-A (200 GeV) 
V 0-Gu, 0-U 
♦ S-U 
\|/' "S-U 



+ 



+ 



-■abs 



: 6.9 + 1.1 mb 



1 



3 4 
L(fm) 



Figure 32: Jjil) and t\}' production in nuclear collisions, compared to pre-resonance ab- 
sorption in nuclear matter |100| , |104| . 



proportional to A^^"^. The transverse energy E'y carried by the secondaries is measured 
experimentally, together with the J/ip^s. We thus have 



Erih) = q 



(73) 



the proportionality constant q has to be determined on the basis of the given experimental 
acceptance. Once it is fixed, we have to check that the collision geometry (the measured 
relation between Et and the number of spectator nucleons, the measured ii^r-distribution) 
is correctly reproduced |5^. When this is assured, we can check if the J/ip production in 
0—Cu, — U and S—U, as measured by the NA38 experiment at CERN over the past ten 
years | ]10CI| , |104| |, shows anything beyond the expected pre-resonance nuclear absorption 



ccg- 



-N 



6.9 ±1.1 mb determined from p — A interactions. 



with the cross section a, 

The answer is clearly negative, as seen in Fig. ^ for the integrated cross sections as well 
as for the centrality dependence of S — U collisions. We consider here the J/ ip survival 
probability, i.e., the measured production rate 5'exp, normalised to the predicted rate 
Sqi, including pre-resonance suppression with the given absorption cross section. To 
compare all data in one figure, we have parametrized the different collision configurations 
in terms of the path length L of the charmonium state in the nuclear environment. As 
a rough estimate, L ^ (3/4) i?^ in p — A collisions, with Ra = 1.13 A^^^ for the nuclear 
radius; the factor (3/4) takes into account the average over centrality. Similar estimates 
can be obtained for A — B collisions; however, the Glauber formalism allows a direct 
determination of this path length at each impact parameter, and we shall use the values 
of L thus determined and then related to through Eq. (|73]). 

All nuclear collisions measured by NA38 thus show only what is now called 'normal' 
J/ip suppression, i.e., the pre-resonance suppression already observed inp — A interactions 
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normalized to pre- resonance absorption in nuclear matter |100, 107 
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Figure 34: The i^^-dependence of J/ip production in Pb — Pb collisions, normalized to 



pre-resonance absorption in nuclear matter |108|. 
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p2 , 101 , 102 |. We thus have to ask if the A — B colhsions studied by NA38 lead to any 



produced secondary medium at all. The behaviour of if)' production (included in Fig. |32D 



shows that this is indeed the case ||100| , |105|| : in S — U collisions, there is ip' suppression 



beyond the expected pre-resonance nuclear absorption. The collisions thus do form a 
secondary medium, and this can distinguish a ip' (which is suppressed) from a, J /ip (which 
remains unaffected). Since the former is suppressed, the latter unaffected, the medium 
does not signal colour deconfinement. 

The observed ip' suppression should thus be accountable in terms of a medium con- 



sisting of hadrons, and such a 'hadronic comover' description ||1U6[| is indeed found to 



work [^. It involves as basic quantities the -i/^'-hadron dissociation cross section and the 



transverse density of hadronic comovers. Since the ip' is very loosely bound (it lies only 
60 MeV below the open charm threshold, see Table 2), it seems reasonable to assume that 
the cross section attains its geometric value of about 10 mb soon after the threshold. The 
density of wounded nucleons is given as function of collision centrality by the mentioned 
Glauber formalism, and the observed rate of about 1.5 secondary hadrons per wounded 
nucleon determines the comover density. The resulting form for ip' suppression is found 
to agree quite well with the existing data |^ . 



5.3 Anomalous J/ip Suppression 

The suppression of J /ip production in nuclear collisions, from p — Ato central S — U 
collisions, is thus understood in terms of 'normal' absorption in nuclear matter: the 
pre-resonance charmonium state is dissociated by collisions with the primary target or 
projectile nucleons. The J/ip as probe of the produced environment indicates "no decon- 
finement". 

This situation came to an end with the advent of truly heavy ion experiments. In 
Pb-Pb collisions at the CERN-SPS, the NA50 Collaboration [^, ^ observed a further 
'anomalous' suppression of about 30% for the overall J/ip production rate (Fig. p3D . 
Looking in more detail, it was observed that peripheral Pb—Pb collisions show only normal 
suppression (Fig. 3^). Then, at a transverse hadronic energy Et — 40 GeV (corresponding 



to an impact parameter b of about 8 fm), the anomalous J/^p suppression sets in quite 
suddenly, with a drop of 20% for an Et increase of less than 20 GeV, or a decrease in b 
of about 1 fm. For very central collisions, E^ > 100 GeV or 6 < 2 fm, something like 
a second drop appeared. The most central collisions thus show a suppression by more 
than a factor two beyond the normal pre-resonance absorption. The general suppression 
pattern, combining all data from p—p to central Pb — Pb collisions, is illustrated in Fig. 

where the J/ip survival probability is shown as function of the number of participant 
nucleons. What conclusions can be drawn from this behaviour? 



5.3a Deconfinement vs. Hadronic Comovers 



Over the past decade, quarkonium suppression mechanisms have become classified 
into two categories. The suppression could be due to colour deconfinement, as originally 
predicted ||6^. To verify this, one has to exclude what are generally called 'conventional' 
mechanisms, covering all possible ways to suppress J/ip production without invoking 
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deconfinement. 



There are a considerable variety of possible conventional suppression 

absorption by produced secondary hadrons 



schemes: absorption in nuclear matter ||111 
(i.e., by the 'hadronic comovers' mentioned above ||106| , |110|| ), nuclear shadowing of parton 



distribution functions ||112|| , medium- induced energy loss of the incident partons leading 
to cc production ||113| ] or more general in-medium modifications of the cc production 
for an extensive review, see [115]. Since all models involve several more 



process 



or less adjustable parameters, the task of confirming one and ruling out all others seems 
quite difficult indeed. It is greatly simplified, however, by two basic qualitative features 
of deconfinement, which distinguish it from all conventional approaches. 

Deconfinement is a critical phenomenon, with a well-defined onset for colour conduc- 
tivity specified by theory. There is a thermodynamic region for confinement, a critical 
point, and then a region of deconfinement. All other approaches are always present in 
smoothly varying degrees. To illustrate the point, we consider two simple models, one for 
hadronic comover absorption, as representative and most intensively studied conventional 
model pO|] 



lis 



the other for colour deconfinement 
In the comover picture, the J/ip is dissociated with a cross section (Xcom in a medium 
of produced hadronic secondaries; the density n{e) of this medium is determined by the 
initial energy density e of the produced environment. The medium subsequently cools 
off and stops affecting the J/ip's when its density drops to a freeze-out value nj. The 
resulting J/ip survival probability becomes [|116|| 



'5}M(e) = exp{-(Tcom^(e)roln[n(e)/n/]}, 



(74) 



where tq ~ 1 fm denotes the formation time of the medium. According to Eq. (|46|) , the 
cross section for J/ip dissociation in the comover medium is expected to be very small, 
ruling out significant contributions from comover absorption. An essential ingredient 
in all comover models is thus the assumption that the short distance QCD arguments 
leading to Eq. (8) become applicable only for much heavier quarks and cannot be used 
for charmonium states.^ The cross section acom is then treated as adjustable parameter, 
as is the freeze-out density Uf. Typical fit values [ 120 lead to (Xcom — 4 - 5 mb. We shall 
return shortly to the latest and most detailed comover fits. 

In contrast to comover absorption, deconfinement sets in at a critical point ec] a simple 



model |[L19|] leads to 



idee 



e(ee - e) + e(e 



ecj 1 - 



9/4 



(75) 



The power of the variable e/ec is determined by the energy density profile in the collision, 
since the hot interior melts J/ip^s, while they survive in the cool rim. The specific form 
( [75| ) arises for nuclei of constant nuclear density [ 121 |. More recent calculations of this 



type ||122|| have led to a very similar suppression pattern, while percolation arguments, as 
will be seen below, lead to a still more abrupt onset of deconfinement suppression. The 

■^The behaviour of a^^ {J / ip — h) for low collision energies can be measured directly either by shooting a 
nuclear beam at a hydrogen target ('inverse kinematics' experiment |117|) or by annihilating antiprotons 



in a nuclear target |118[ . Hence this question can be settled empirically. 
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abrupt onset of J / ip suppression in Eq. (|75|) shows a fundamental qualitative difference to 
the smooth monotonic decrease of the comover suppression pattern Eq. ([74D; an explicit 
illustration will be given shortly. It should also be noted that hadronic comover absorption 
would have to be present already in S — U collisions, thus ruling out pure pre-resonance 
absorption here. The behaviour seen in Fig. p2| seems difficult to reconcile with this. 
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Figure 35: J/ip production in A~ B and Pb — Pb collisions, normalized to pre-resonance 



absorption in nuclear matter, as function of the number of participating nucleons ||108| , 



So far, we have addressed the fate of a well-defined J/ip in a. strongly interacting en- 
vironment. However, sufficiently energetic nuclear collisions produce the entire spectrum 
of charmonium states; below the open charm threshold, they produce r]c{lS), J/ip{lS), 
Xco(l-P)) Xcii^P), Xc2(l-P) and ip'{2S). All higher excited states have significant decay 
rates into the J/ip, and so the observed J/ip's can be either directly produced, or they 
can come from the decay of higher excitations. Since the decay times of these excitations 
are very long (more than 1000 fm for the Xc), the medium affects the intermediate higher 
excitations, not the final decay product J/ip. 

In nucleon-nucleon collisions, about 60% of the observed J/ip's are produced directly, 
about 32% come from Xc decays, the remaining 8% from ip' decays. To simplify matters, 
we neglect the ip' component for the moment and assume 2/3 direct production and 1/3 Xc 
decay. Both colour screening and gluon dissociation indicate that the Xc will be destroyed 
at a lower energy density than the direct J/ip's. As second basic qualitative feature, colour 
screening thus predicts sequential J/ip suppression, with a multi-step survival probability 



54| , |121|| . For J/ijj suppression, the simple analytic form (|75D now becomes 



odec I 



9(6^ - 6) + 9(6 




9(e^-e) + 9(e-e 




(76) 
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where ~ ec and > ec denote the thresholds for x and direct i}) dissociation, re- 
spectively. This form of the resulting functional behaviour of deconfinement suppression 
(Fig. |3^) was given eight years ago [|119|| ; it clearly illustrates the quantitative difference 
to suppression by comover absorption. Today we find it to be remarkably similar to the 
suppression pattern now observed in Ph — Ph interactions (Figs. 03 and BSf). 



The first Ph — Ph results provided by NA50 [123] were of much lower precision than 
those now available. This stimulated various comover approaches to attempt a unified 
description of p — A, S — U and Ph — Ph data. The monotonic onset of any comover 
suppression requires the presence of some anomalous suppression also for S — U, and it 
excludes any thresholds. Hence even the first, relatively pronounced onset of suppression 
in the Ph — Ph data posed problems. Nevertheless, qualitatively the comover pattern also 
leads to a decrease with increasing energy density, and so more precise data were needed 



for any reliable conclusion. Sucessive NA50 runs have provided this ||107| , |108|| , and the 
given experimental situation today, with a clear two-step suppression pattern, excludes 
all conventional scenarios. In Fig. we show as illustration the results of one of the most 



general comover approaches, based on the dual parton model | |125| |; it includes comover 



production dependent on the number of collisions as well as on the number of participant 
nucleons, and it also takes into account J/ tp'-J/ transitions. In spite of the resulting 
number of adjustable parameters, it is clearly incompatible with the latest data. All other 
conventional models share this fate: none can account for all the available data, including 
S — U and Ph — Ph collisions; in particular, none can produce a mechanism leading to any 
kind of non-monotonic step structure. This provides the basis for the conclusion reached 
by the NA50 collaboration [|108|| : The onset of the anomalous J/ip suppression is the first 
clear observation of a threshold effect in heavy ion collisions and can be considered as a 
strong indication of the production of a deconfined quark-gluon plasma phase in central 
Pb-Pb collisions. 

5.3b Towards an Understanding of Deconfinement 

The observed J/ip suppression thus rules out conventional explanations; moreover, it 
is qualitatively in accord with the behaviour expected from colour deconfinement. Nev- 
ertheless, the details of the onset of deconfinement are not yet clear, and in particular 
the crucial collision variable triggering this onset is not yet unambiguously identified. We 
present in this section one possible approach, based on parton percolation |^ p.28| . 



In a high energy nuclear collision viewed in the overall center of mass, the two Lorentz- 
contracted nuclei quickly pass through each other. After about 1 fm at the SPS and only 
about 0.1 fm at RHIC, they have separated and left behind a partonic medium as a 
potential QGP candidate. Inside this medium, primary nucleon-nucleon collisions have 
left charmonium states as probes, whose fate can provide information about the nature 
of the partonic medium. 

The average number of produced partons will depend on the number of nucleon- 
nucleon collisions and/or the number of participating nucleons. Up to SPS energies, 
the number of wounded nucleons appears to be the main determining factor. This is 
understandable in terms of the dilated and hence large soft parton formation time as seen 
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Figure 36: J/ip suppression by deconfinement compared to that by hadronic comover 
absorption ||119|| . 
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Figure 37: J /if) production in Ph — Ph collisions []6T| , |124|| , compared to the dual parton 
model comover approach 
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in the rest frame of either nucleus. Interference and cancellation effects of the Landau- 
Pomeranchuk type thus prevent soft parton emission at each nucleon-nucleon collision. At 
higher energies, large additional contributions can come from hard partons (minijets or 
jets) with short formation time, and these will be proportional to the number of collisions. 
Even at SPS energy, nucleon stopping and secondary multiplicities at mid-rapidity seem to 
increase somewhat with increasing mass number A, and this could be the onset of collision 
dependent effects. The s-dependence of the mid-rapidity multiplicity is presumably a 
combination of the s-dependence of the parton distribution function and the onset of 
significant hard hadron production. We shall try to include these dependences here in 
rather phenomenological terms, without specifying their origin. 

The number of gluons of transverse momentum kx emitted by a wounded nucleon 
can be obtained from the nucleonic gluon distribution function g{x) from deep inelastic 
scattering studies, 

/ dNg 



dydkj 



^ ^ xg{x)m) (77) 

where x ~ kT/y/s. Here ^/s is the incident collision energy, and f{k^) is normalized to 
unity. Using the MRS-H form of the gluon distribution function ||126|] and integrating the 
resulting Eq. ( [TTD over kx then gives us {dN/dy)y=Q ~ 2 for SPS and 4 for RHIC energy. 

Each gluon of transverse momentum kx has an effective transverse size r ^ k^^, and 
nucleus- nucleus collisions provide many such gluons overlapping in the transverse plane. 
We want to study their clustering behaviour. To simplify matters, we shall consider 'av- 
erage' gluons (rather than averaging results with f{k^)). From the transverse momentum 
dependence of Drell-Yan dilepton production through quark-antiquark annihilation, or 
from that of J/ip production through gluon fusion, one finds that the effective intrinsic 
transverse momentum of gluons is (kr) — 0.75 GeV/c, leading to an average transverse 
parton radius r ~ 0.27 fm. We now want to consider the deconfinement pattern for such 
systems. 

The collisions distribute the gluons over the transverse area according to the nucleon 
distribution within the nuclei. We thus have to distribute discs of radius r ~ 0.27 fm in the 
transverse plane and determine the overall nucleon density for which colour conductivity 
sets in, i.e., for which there appear clusters which reach the dimension of the entire 
system. Hence we again encounter percolation as crucial feature for deconfinement; here 
it occurs as a two-dimensional phenomenon, since djTiamics aligns the gluon motion along 
the collision axis. If we denote with ric the critical gluon density for the production of 
a deconfined medium, then as in section 2.4, it appears natural to identify the onset of 
quark-gluon plasma formation with the percolation point, which is p5| 



ric = ^, (7^ 



where vrr^ is the transverse gluon area. If we consider gluons of transverse momentum 
(kx) = 0.75 GeV/c as characteristic for nucleon-nucleon collisions at present CERN-SPS 
energies, then the threshold obtained from Eq. (|78D is 7ic — 5 fm~^. The deconfined 
medium thus obtained consists of gluons freed of their confinement constraints, but it 
does not yet need to be in thermal equilibrium. 
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Figure 38: Average cluster density nci(n) (left) and average fractional cluster size 
Sc\{n)/Stota.\ (right) as function of the overall density n of discs, for r/R = 1/20; in 
(b), the derivative of Sci{n) / Stotai with respect to n is also shown (dotted line). The 
percolation point in the limit r/R —>■ is indicated by Up. 



As a prelude, we assume N discs of radius r to be randomly distributed on a flat 
surface of radius R » r and study the average cluster density nci(n) and the average 
cluster size S'ci(n) as a function of the overall density n = N/nR^. The result is shown in 
Fig. |38|, where we have normalized all quantities to the inverse disc area l/vrr^ in order 
to make them dimensionless. We note that the derivative of the cluster density peaks 
sharply at a certain density; in the 'thermodynamic' limit R —>■ oo, this percolation point 
is known to be Up = 1.175 . We see in particular that if we want to reach a certain cluster 
density n^, attained at a certain overall n, then at this point the average cluster S'ci(n) 
also has a certain finite size. In other words, requiring a specific density for the onset of 
a new phase implies that this onset occurs for a specific, finite size of the system. The 
percolation point, for example, is reached when about 50% of the surface is covered by 
discs. 

Turning now to nuclear collisions, we have to distribute the partonic discs not on a 
fiat surface, but instead in the transverse plane according to the nucleon distribution 



determined by the profiles of the colliding nuclei ||127|| . Moreover, we have to allow 
collisions at different impact parameters, and we want to consider different A—B collisions. 

In Fig. |3^, the cluster density Ud is shown as function of the overall density riparton 
of 'average' partons of radius r = 0.27 fm, for different centralities of various A — B 
combinations at SPS energy, with 2.1 partons per wounded nucleon. For each A — B, the 
point at highest riparton corresponds to central (impact parameter 6 = 0) collisions, the 
one at lowest nparton to the most peripheral collisions. We note that by varying A — B 
and centrality we produce an essentially universal cluster density nci(nparton)- Now we 
assume that deconfinement occurs at the percolation point, which can be determined by 
studying the average fractional cluster size 5'ci/S'totai- In Fig- ^0], we see the result for 
Pb — Pb collisions. At the density for which the derivative peaks, nparton — 4.2 fm~^, the 
cluster density reaches its critical value rid — 6 fm~^. The same critical value is found for 
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Figure 39: Cluster density vs. parton den- 
sity for different centralities and different 
A — B configurations at SPS energy. 



Figure 40: Fractional cluster size vs. par- 
ton density, together with its derivative, 
for Pb — Pb collisions at SPS energy. 



Sn — Sn and U — U collisions at SPS energy, even though this value is attained at quite 
different centralities in the different configurations {b ~ 4.5 fm for Sn, 8.5 fm for Pb and 
11.5 fm for U collisions). 

In Fig's. ^ and we show the result of this procedure for RHIC energy, where we 
have 4.0 gluons per wounded nucleon. Again the critical cluster density at deconfinement 
is found to be rid — 6 fm~^. Here we have to consider Cu~Cu collisions in order to 
study the onset of deconfinement, since for Pb — Pb all possible centralities are above the 
threshold. 

We thus conclude that the critical cluster density is, as expected, a universal quantity, 
independent of the choice of A—B, of centrality, and of incident energy. The fractional size 
of the deconfined cluster at threshold does, however, depend on the collision configuration: 
for Pb-Pb at SPS, we have Sd/Stotai ^ 0.45, while Cu-Cu at RHIC gives 0.70. In other 
words, at the higher RHIC energy, deconfinement sets in for a larger bubble than at the 
lower SPS energy. 

We further note that at SPS energy, S — S and S — U are below the deconfinement 
threshold even for the most central collisions, while at RHIC even the most peripheral 
Pb—Pb collisions are above the threshold. To study the threshold, we thus require heavy 
nuclei (such as Pb) at the SPS, lighter nuclei (such as Cu) or lower incident energy at 
RHIC. 

Having specified deconfinement through percolation, we now turn to J/ ip suppression 
as a probe for its onset. As noted above, the Xc is dissociated essentially at deconfinement, 
so when the cluster density reaches n™* ~ 6 fm~^, all Xc's which are inside the percolating 
cluster disappear. The amount of Xc suppression thus depends on the fractional cluster 
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Figure 41: Cluster density vs. parton den- 
sity for different centralities and different 
A — B configurations at RHIC energy. 
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Figure 42: Fractional cluster size vs. par- 
ton density, together with its derivative, 
for Cu—Cu coUisions at RHIC energy. 



CO 



1.2 
1 

0.8 
0.6 
0.4 
0.2 







1 1 1 


- • 1 

b = 1 1 fm 




b = 


SPS : Pb-Pb 


1 1 



'parton 



Figure 43: The J /if) survival probability as function of the parton density for Pb — Pb 
collisions at SPS energy; the dotted line corresponds to a medium of uniform parton 
density, the solid line to collisions with parton densities determined by the profiles of the 
colliding nuclei. 
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Figure 44: The J/ip survival probability as function of the parton density for Ph — Ph 
(left) and for Cu — Cu (right) collisions at RHIC energy; the dotted line corresponds to 
a medium of uniform parton density, the solid line to collisions with parton densities 
determined by the profiles of the colliding nuclei. 



size at deconfinement, which is about 45 % for Ph — Pb at the SPS, about 70 % for 
Cu - Cu at RHIC and effectively 100 % for Ph - Ph at RHIC. At deconfinement, the 
corresponding fractions of J/ ip production through intermediate Xc states should thus be 
suppressed. 

As already noted, the dissociation of directly produced J/ip^s requires a a larger screen- 
ing mass and thus higher density than available at the point of deconfinement. On a 
microscopic level, J/tp break-up becomes possible only for harder than average gluons at 
deconfinement. In principle, the relevant quantities can be determined by lattice QCD 
studies; however, this requires computer capabilities which are being reached only now. 
For the moment, we therefore choose v^^^ = 0.22 fm and n'^l^ = 7.8 fm~^; this corre- 
sponds to a ratio of transverse dissociation energy densities eri^J I ip) / ^TiXc) — 1-6, in 
accord in studies combining potential theory and lattice results ||6^, |121|| . With these 
parameters, one obtains 20% direct J/ip suppression in Ph—Ph collisions at SPS energy, 
once the required cluster density is reached. The corresponding fractions for RHIC are 
30% for Cu-Cu and 35% for Ph-Ph collisions. 

Schematically, one thus obtains a suppression pattern as illustrated in Fig. where 
it is assumed that 40% of the observed J/ip's come from Xc decay and 60% are produced 
directly. For clarity purposes, we ignore here the small fraction (about 8%) coming from 
ip' decay. The dotted curve corresponds to the suppression which would take place in 
a uniform medium of precisely specified parton density: at the deconfinement point, all 
Xc are dissociated, so that the corresponding fraction of decay J/ip's is gone; when the 
density for direct J/ip melting is reached, these disappear as well, leading to complete 
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J/ip suppression. In actual nuclear collisions, the medium is not uniform, with denser 
'inner' and less dense 'outer' regions in the transverse plane. Suppression now occurs 
only in the fraction S'^^/S'totai for the Xc part and Sj/^/ Stotai for the direct J/ip part, 
where S^^ and Sj/^ denote the clusters in which the respective dissociation density is 
reached. Since these fractions depend on energy and nuclear geometry, the suppression 



curves are different for different experimental configurations. In Fig. 43, the result is 



shown for Pb — Pb collisions at the SPS, from impact parameter 6 = 11 fm to 6 = fm. 
The difference between this curve and the one for a uniform medium thus reflects the 
surviving J/ip^s produced in the less dense outer regions. 

In Fig. ^ similar calculations are shown for Cu — Cm and Pb — Pb at RHIC energy. 
Here we note in particular that Pb—Pb collisions are for all meaningful centralities {b < 11 
fm) above both the Xc and the direct J/ip threshold, so that we get a smooth anomalous 
suppression increasing from about 60% at 6 = 11 fm to about 90% at 6 = 0. Combining 
this suppression with the 'normal' pre-resonance absorption in nuclear matter, we thus 
predict for central Pb — Pb collisions at RHIC a. J ftp production rate of less than 5% the 
corresponding unsuppressed rate (excluding possible B decay contributions). It should 
also be noted that this result is based on twice the number of gluons per wounded nucleon 
at RHIC, compared to the SPS value. A larger increase, based on a possible larger hadron 
multiplicity at RHIC, would lead to more J/ip suppression. Similarly, we assume an 
average number {dNg/dy)y=o of gluons per wounded nucleon. It is conceivable that very 
central collisions reach into the tail part of the multiplicity distribution, with a larger 
number of hadrons and hence also gluons. This would lead to a larger suppression for 
very central collisions. At the SPS, the NA50 collaboration can check if the basis for 
this exists. By combining measurements of hadron multiplicity, transverse energy Ej- and 
forward energy Ezbc, it is possible to study the number of hadrons per wounded nucleon 
as function of Et and check if there is an increase at highest Et values. 

Finally we can then turn to a comparison of our results on parton percolation to 
the i?T-dependence of the actual data, which contain an additional smearing due to the 
fact that a given E^ bin corresponds to a range of impact parameters and hence parton 
densities. Including this effect in the standard way [^, we obtain the result shown in 



Fig. ^ Included in this figure are the 1997 and 1998 data ||124| , |6T| with minimum bias 
determined Drell-Yan reference. While reproducing the overall trend, our results clearly 
show deviations from the data in detail. As mentioned, the levelling-off of our curve at 
high Et would have to be modified if an increase in the multiplicity per wounded nucleon 
at high Et should be observed. Another possible modification could enter through a 
density-dependent charmonium dissociation. We have here assumed that all Xc's melt 
once the critical density is reached. Allowing a partial survival chance, which decreases 
with increasing density, would lead to a steeper drop of the suppression with nparton as 
well as with Et- One possible source for such an effect could be the finite life-time of the 
deconfining medium in actual nuclear collisions. 

Summarizing this section, we note that parton percolation provides a consistent frame- 
work to study the onset of deconfinement as well as the onset of charmonium suppression 
as a signature of the transition. It is clear that the compatibility of the parton basis with 
other, soft hadron production processes has to be checked; we return to this in section 6. 
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Figure 45: J /ip production in Ph — Pb collisions []6T1 , |124|] , compared to colour deconfine- 
ment in the parton percolation model [12S]. 



5.4 Pt Dependence 



In this section, we want to consider the dependence of anomalous J/ip suppression 
on the transverse momentum of the J/ip. It is well known that the transverse momenta 
of secondaries from hadron-nucleus collisions quite generally show a pT-broadening. This 
means that the universal parameter A ~ 0.15 GeV in Eq. ( ppj ) now depends on the size 
of the target nucleus. For secondary hadrons, this is the Cronin effect [ |129|| ; a similar 



behaviour is observed also in Drell-Yan and charmonium production ||9^. The natural 
basis for all such broadening is initial state parton scattering, and it was in fact shown 
some time ago ||130|| - ||133|| , [p.l9|| that this describes quite well the pT-dependence observed 
in J/ip production from p—A to central S — U collisions. 

Consider J/ ip production in p — A collisions, assuming gluon fusion as the dominant 
process for the creation of a cc pair. Parametrizing the intrinsic transverse momentum 
distribution /(gr) of a gluon in a nucleon as 
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we obtain by convolution for the transverse momentum distribution Fpa{Pt) of the re- 
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with 

{Pt)pA = {qDa + {ql)p. (81) 

The quantity 

8pA = {Pt)pA — {Pt)pp = {Qt)a — {Qt)p (82) 

is thus a suitable measure for the observed nuclear broadening. Note that the Gaussian 
form of Eq. (|8^ is obtained from Eq. (^) for P| << rrij^^. For J/ip production in pp 
collisions, with the universal ~ 0.15, this leads to (P^) ^ 2mj/^/X ~ 1 GeV^, which 



is slightly lower than the observed value of 1.23± 0.05 ||97 |. 

Assume now that in the passage of the projectile proton through the nuclear target, 
successive interactions broaden the intrinsic momentum distribution of the corresponding 
projectile gluon which will eventually fuse with a target gluon to form a J/ ip. If the process 
of Pt broadening during the passage is a random walk, then the relevant parameter of 
the Gaussian distribution (^) becomes 

SpA = N%, (83) 

where is the average number of collisions the projectile undergoes on its passage 
through the target up to the fusion point, and 6o the average broadening of the intrinsic 
gluon distribution per collision. 

In nucleus-nucleus collisions, a corresponding broadening occurs for both target and 
projectile gluon distributions; here, however, measurements at fixed transverse hadronic 
energy Et can determine the broadening for collisions at a given centrality. Hence at 
fixed impact parameter b we have 

SABib) = {P^ABib) - (P|)pp = N^'^ib) 5o, (84) 

with N^^ib) denoting the average number of collisions for projectile nucleons in the target 
and vice versa, at fixed b. N^^{b) has a maximum at small b and then decreases with 
increasing 6; for a hard sphere nuclear model, it would vanish when b = Ra + Rb- 

In Glauber theory, the quantity N^^{b)/a can be calculated parameter-free from the 
established nuclear distributions |127|; here a denotes the cross section for the interaction 



of the nucleon on its passage through the target. We shall determine (t5o from data, so 
that (J never enters explicitly. Once cr^o is fixed, the broadening by initial state parton 
scattering is given for aS[ p — A and A — B interactions. For Drell-Yan production (with 
quarks instead of gluons in the partonic interaction), this would be the observed effect, 
since the final state virtual photon does not undergo any further (strong) interactions. 
A produced nascent J / ip will, however, experience pre-resonance nuclear absorption; this 
suppresses J/ip^s produced early along the path of the projectile, since they traverse more 
nuclear matter and hence are absorbed more than those produced later. As a net result, 
this shifts the effective production point to a later stage. In p — A collisions, a Drell-Yan 
pair will on the average be produced in the center of the target. In contrast, nuclear 
absorption shifts the average cc production point further down-stream. This effectively 
lengthens the path for initial state parton scattering and hence increases the resulting 
broadening. 
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The transverse momentum behaviour of normal J/ip production in nuclear collisions 
is thus a combination of initial state parton scattering before the production of the basic 
cc state, and pre-resonance nuclear absorption afterwards; both lead to a broadening of a 
Pr-distribution. A further broadening could come from elastic random walk scattering of 
the charmonium state itself in nuclear matter; however, the effect of this can be taken into 
account by fitting cr^o to the data. The essential task is then to calculate the number of 
collisions per cross section, N^./ cr, for p —A and A—B interactions, taking into account the 
effect of pre-resonance nuclear absorption. The result contains as only open parameter 
the quantity a5] for details, see ||136 |. 

Ideally, one would use p — A data to fix a6o; the broadening for A — B interactions 
would then be fully predicted. Unfortunately there are p —A data only for three values of 
A [p7| , p9| , and these have rather large errors. It is therefore better to attempt a consistent 



description of all existing p — A |134|| and A — B data |p.34|| , up to central S — U, in 
terms of a common ct^q. Using {Pt)pp = 1-23 ± 0.05 GeV^, a best fit to p — A, — Cu, 
0-U and S-U data [|l||] gives a5o = 9.4 ±0.7, with a x^/d-f. of 1.1. In Fig. || we show 
the resulting description of the S — U data. 

We then turn to Pb—Pb collisions; the corresponding "normal" transverse momentum 
behaviour is shown in Fig. ^7\. Also shown here is the amount of broadening obtained from 
initial state parton scattering alone; we see that nuclear absorption increases {P^)pbPb 
by about 15 %. Since the functional forms in the two cases are quite similar, nuclear 
absorption can be simulated by choosing a somewhat larger 5q. The basic feature of 
normal absorption thus remains the monotonic increase of (P^) with Et-, even though the 



collision geometry makes this slightly weaker for Pb—Pb than for S—U interactions ||135 



The onset of anomalous suppression results in a striking modification of this pattern. 
If the J/ip^s in the hot interior of the medium produced in the collision are suppressed. 
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Ej (GeV) 

Figure 47: Normal Pt broadening in Ph — Pb collisions with (solid line) and without 
(dashed line) pre-resonance nuclear absorption. 




Figure 48: Pt broadening in Pb — Pb collisions. The solid line shows normal broadening, 
the dotted line the effect of anomalous Xc ^-nd t/j' suppression, and the dashed line the 
effect the onset of direct J/ip suppression. 
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then this will reduce their contribution from the part of phase space leading to the most 
broadening. To illustrate the effect, we assume suppression by deconfinement; in this case 
the result is readily calculable [p.l9| , p.32|| . To be specific, we assume deconfinement to 
start once the number of collisions per wounded nucleon exceeds a certain critical value 

N^{b,s) 

This defines for each collision configuration cool outer and hot inner regions, resulting in 
charmonium survival or dissociation, respectively. In Fig. 4S we compare the resulting 
patterns for several value of Kc to that of normal absorption; the change of behaviour 
is qualitatively evident and due solely to anomalous J/ip suppression. Moreover, abrupt 
onsets of anomalous Xc and direct J/ ip suppression would also show up here; they have not 
yet been included in Fig. E8|. The observation of such anomalous transverse momentum 



behaviour for J/^ suppression would thus constitute another tool to study the onset of 
deconfinement in Pb — Pb collisions. 



6. SOFT HADRONIC PROBES 

In this final section, we want to address the possible supporting evidence which can 
be obtained by measuring soft hadronic probes, i.e., mesons and baryons made up of 
the light u, d and s quarks and produced at low transverse momenta. Since there is no 
large intrinsic mass or momentum scale, their production occurs over comparatively large 
space-time regions determined by Aqcd and cannot be described perturbatively. The 
large intrinsic size of soft hadrons requires relatively low energy densities to make their 
existence conceptually meaningful: five mesons need a volume of at least roughly five 



meson volumes in order to exist ||53[. As discussed in section 4.1, this implies that they 
are formed late in the collision evolution, when the high initial energy density has dropped 
to a sufficiently low value eu < 0.3 — 0.5 GeV/fm^, taking the energy density of a typical 
hadron as indicator. 

What features of the produced medium can one then study with the help of soft 
hadronic probes? It is clear, as also noted in section 4.1, that they provide direct infor- 
mation about the state of the medium at the energy density eh- In addition, one may 
hope to obtain 'retrospectively' also information about earlier stages at higher e. In the 
following, we shall first consider two specific instances of such retrospective probing and 
then discuss the study of the medium at hadronization. 

6.1 Retrospective Probes 

The onset of anomalous J/ip suppression occurs at a certain associated transverse 
energy of secondary hadrons. We would like to translate this into a value of the corre- 
sponding primordial energy density. In principle, this can be done in the same way as 
was used to derive Eq. (|36D in section 3.2. It is clear, however, that the result will be 
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model dependent. Instead of the assumed free flow leading to Eq. (0), one could assume 
an isentropic flow ||137|| , leading to 

as initial entropy density. It turns out that with Eq. (|35|) for the multiplicity and for an 



ideal gas of massless partons, the two different expansion patterns result in very similar 
initial energy densities up to ~ 200 GeV. At much higher y^, isentropic expansion 
leads to a somewhat higher energy density, but even at LHC energy, the difference is only 
about 20 % 

To obtain an estimate of the energy densities in the region of anomalous J/?/^ sup- 
pression, we thus make use of Eq. (PU]). For the 5% most central Ph — Ph collisions, one 
finds on the average about 800 secondaries at central rapidity; their average energy is 
about 0.5 GeV ||144|| . This leads to an initial energy density of about 2.8 GeV/fm^. Note 
that this value corresponds to e averaged over the collision profile; the 'hot' central value, 
associated to the region of highest density of participating nucleons, will be around 30% 
higher, around 3.7 GeV/fm^ |52]. This is then the value of e at the point of maximum 



J/V' suppression in the SPS data. To study in detail the onset of deconfinement, the mul- 
tiplicity of soft secondaries must be known as function of centrality. Until a detailed study 
of this is available, we can estimate the result by assuming the e(6) (and hence e{ET)) 
is approximately linear in Et- This leads to ec — 2 GeV/fm^ at the onset of anomalous 
J/ip suppression, which presumably occurs at the onset of Xc dissociation and hence at 
deconfinement. Such a critical energy density is in reasonable agreement with the values 
predicted in section 2. 
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A second and particularly interesting retrospective probe could well be the strangeness 
abundance in nuclear collisions [|139|| . We had seen in section 3.1 that the relative abun- 
dances of the secondary species in p — p collisions were determined by the freeze-out 
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temperature of a hadronic resonance gas, with the restriction that strangeness pro- 
duction was suppressed by a factor 7^. In nuclear colhsions, the non- vanishing baryon 
number density of the system requires that a corresponding analysis specifies temperature 
Th and baryochemical potential /i^ at freeze-out, and in addition a possible strangeness 
suppression factor 7^. In central nuclear collisions, strangeness production is found to be 
much more abundant than inp —p collisions at the same ^/s. As illustration, we show in 
Fig. HSIthe K/ti ratio in Ph — Ph collisions at the SPS as function of the number of partic- 
ipating nucleons [|141| , |142|| . If we assume, as discussed in section 3.1, that the suppression 
factor 7s is determined by the relative formation rate of the heavier strange quarks in an 
earlier partonic medium of higher temperature T > T^, then 7s (T) could provide a ther- 
mometer for this early temperature T. The behaviour seen in Fig. ^ would thus indicate 
a significant increase with centrality for the initial temperature of the partonic medium 
in Pb — Pb collisions at ^/s = 20 GeV. Comparing Fig. ^ with Fig. 0, we see that the 
strangeness increase with centrality in nuclear collisions is very similar to that found as 
function of ^/s or as function of multiplicity at fixed y/s in the p — p/p — p interactions. 
This supports an increase of the initial energy density as the underlying mechanism for 
increased strangeness production; but further comparative studies seem necessary to fully 
specify the origin of the effect. If the model for 7s (T) discussed here is correct, we would 
expect that with increasing collision energy, at RHIC and LHC, 7s — 1. We also note 
that in accord with the interpretation of 7s as thermometer, there is so far no indication 
of a threshold or sudden increase in the relative strangeness production. 
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6.2 Hadronisation and Freeze-Out 



After strongly interacting matter has cooled off to reach the hadronization point, 
its freeze-out into free hadrons should be a universal phenomenon, depending only on 
the freeze-out temperature and baryon density. In view of the expansion of the produced 
medium, this universality can only be 'local', e.g., volume elements moving relative to each 
other can lead to different effective temperatures. In general terms, such an expansion 
scenario leads intuitively to three specific predictions ||143| 

• relative hadron abundances specified by a equilibrated hadronic resonance gas, most 
likely with 7^ ^ 1, 

• a mass-dependent transverse momentum broadening of hadron spectra due to radial 
hydrodynamic flow, 

• a source size for hadron emission which increases with increased expansion due to 
increased initial energy density. 

All three effects are indeed observed. 

For the relative abundances of the soft hadron species, there are indications for the 
validity of a resonance gas distribution in A — A collisions |[144|| - ||146|| . In Fig. this 
is illustrated for Pb — Pb collisions at the SPS; we note that now 7^ is indeed quite close 
to unity. The latest analysis shows some open questions concerning the universality of 



the resulting values of Th, and 7^ | P-4q |. We note here that all studies so far have 
assumed that for a given collision configuration, such a 'hadrochemical' equilibrium exists 
globally over the entire volume of produced matter. In view of the collision profile of 
participating nuclei, howver, it seems quite possible that such an equilibrium exists only 
locally, making the triplet of values T/j, and 7^ dependent on the distance r from the 
collision axis. Since most expansion models start from such local equilibrium, it seems 
natural that species abundances would also be affected by it. 

In Fig. ^ we see that there is indeed a significant (and mass-dependent) broadening 
in the transverse mass distribution of hadronic secondaries from Pb — Pb collisions at 
the SPS | [I44| |. Clearly it would be interesting to study the kinematic counterpart of a 
resonance gas freeze-out and check if the momentum distributions of secondaries become 
thermal and if the expanding system leads to the expected hydrodynamic flow ||143 |. 



Such an endeavor is confronted, however, by a speciflc problem. Besides collective final 
state effects, there are also initial state effects such as multiple state parton scattering 
which can and will modify the transverse momentum distribution Eq. ( p6D of secondaries 
Hence it is not so easy to separate the possible origins of the observed transverse 
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momentum broadening in nuclear collisions | |144|| and to associate them to specific aspects 
of the produced state. Nevertheless, considerable effort has been devoted to this aspect 
[ |143| . For the purpose of our primary question, colour deconfinement, we note it does not 
seem possible to identify the initial state of the produced medium from the hydrodynamic 



expansion pattern alone [p.47|| . 

Because of expansion and momentum-dependent effects, the source size situation is 
somewhat more complex, although an increase with incident collision energy has been 
observed 111 4811 ; for recent discussions, see [1143, 11441 . 
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All three phenomena listed above are present already in p — p/p — p collisions for either 
increasing collision energy or increasing hadron multiplicity. For hadron abundances, this 
was shown in Fig. O; transverse momentum broadening in hadronic collisions is illus- 



trated in Fig. ^ and the the growth of the source size with multiplicity is shown in ||140 |. 
Evidently this supports strongly the idea that increased initial energy density is the com- 
mon origin of the three effects. However, to identify the different underlying mechanisms 
quantitatively, it seems necessary carry out a comprehensive systematic comparison of 
p — p,p — A and A — A collisions, which so far is lacking. As a result, a full understanding 
of the properties of strongly interacting matter at freeze-out and their origin seems to re- 
quire considerable further experimental as well as theoretical studies. Microscopic models 
based on partonic interactions may well provide a suitable tool to formulate a consistent 
picture. 
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Figure 51: The average transverse momentum of different hadronic secondaries in Pb — Pb 



interactions at the SPS ||144|| . 
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6.3 Parton Cascade Models 



For systems evolving in thermal equilibrium, later stages carry no information of the 
early history. In case of a non-equilibrium evolution, such information can be passed on to 
a certain extent. A study of the onset of deconfinement in a microscopic, non-equilibrium 
space-time picture would therefore be of great help in understanding if and to what extent 
the formation and the properties of the quark-gluon plasma are related to experiments 
observing the different evolution stages. The aim of partonic cascade models |50, 51 1 is to 
describe through perturbative partonic interactions embedded in a relativistic transport 
theory the evolution of a high energy nucleus-nucleus collision. Hence such models seem 
to be a good starting point for a dynamic study of deconfinement and the subsequent 
evolution of the produced systems. 

We had seen above in section 5.3 that the initial state of two colliding nucleons can 
be pictured as colliding beams of confined partons. The incident partons have a distribu- 
tion in intrinsic transverse momentum, leading to some average value (fcr), which in turn 
defines an average transverse parton size tt = l/{kT)- When two nuclei collide at high 
energy, the transverse parton density increases considerably, since now we have a super- 
position of nucleon-nucleon collisions. In this medium, colour screening will destroy the 
association of partons to particular hadrons, since for sufficiently high density of colour 
charges, the colour screening radius becomes much smaller than the typical hadronic 
scale. Hence we expect the onset of deconfinement for some characteristic density or for 
an equivalent screening scale. In perturbative QCD, the method to calculate the colour 
screening mass /i (the inverse of the screening radius ru) is well-known; an extension to 
the non-equilibrium medium provided by cascade models has also been proposed ||149| 
and allows a calculation of the time evolution of the screening mass fi{t) [|151 ]. 

Given fi{t), we have to specify for what value deconfinement sets in; parton cascade 
models themselves does not identify such a point. They just introduce the primary col- 
lisions between the 'confined' partons of the incident nucleons, followed by successive 
interactions between primary as well as produced partons and eventually by the hadroni- 
sation of partons according to some given scheme. There are two possible ways to specify 
deconfinement in such a scheme. One can obtain at any given time the transverse density 
profile of the parton distribution, and a natural onset of deconfinement is the percolation 
point of partonic discs in the transverse plane, as discussed in section 5.3b. A more global 
approach is to calculate the screening mass fj, and identify as critical point fic the value 
obtained from lattice studies. These provide the temperature dependence of /i and in par- 
ticular also its value at the deconfinement point Tc, where one has fic = fJ^iTc) ~ 0.4 — 0.6 
GeV It is not obvious that such an equilibrium value can really be used in the 

non-equilibrium situation provided by the parton cascade model. A good check would be 
to compare in this model the transverse parton density at fic to the percolation value. 

We now briefiy recall the most important features of the parton cascade model; to be 



specific, we consider the version proposed in |50 



The initial nucleus-nucleus system is treated as two colliding clouds of partons, 
whose distribution is fixed by the nucleonic parton distribution functions determined 
in deep inelastic lepton-nucleon scattering, and by the nucleon density distribution 
in the nuclei. 
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• The parton cascade development starts when the initial parton clouds inter-penetrate 
and traces their space-time development due to interactions. The model includes 
multiple elastic and inelastic interactions described as sequences of 2 — 2 scat- 
terings, 1^2 emissions, and 2 — »• 1 fusions. It moreover explicitly accounts 
for the individual time scale of each parton-parton collision, the formation time of 
the parton radiation, the effective suppression of radiation from virtual partons due 
to an enhanced absorption probability by others in regions of dense phase space 
occupation, and the effect of soft gluon interference in low energy gluon emission. 

• Finally, the hadronization in terms of a parton coalescence to colour neutral clusters 
is described as a local statistical process that depends on the spatial separation and 
colour of the nearest-neighbour partons ||150|| . These pre-hadronic clusters then 
decay to form hadrons. 

This model can now be used to calculate characteristic features at different stages of 
the evolution process. As an illustration we show in Fig. |S3 the evolution of the colour 
screening mass as function of the proper time r of the system. It is seen that for sufficiently 
early times, the hot interior of a central Pb — Pb collision leads to h{t) > Hc and thus to 
deconfinement |[151|| . One can now follow the evolution of the model and study features of 
the system beyond the confinement point, such as in- medium hadron mass modifications 
or momentum distributions of soft secondaries ||152|. 




Figure 53: The time dependence of the colour screening mass in central Pb — Pb collisions 



at SPS energy in different radial collision zones ||151 



It is clear that such a description is necessarily model-dependent. Nevertheless, a 
coherent picture of this kind, accompanied by emprirical cross-checks for the different 
evolution stages, provides an excellent visualisation of the expansion process and could 
thus lead to further predictions. It also allows an investigation of problems, such as the 
comparison to p — p and p — A collisions or the onset of thermalisation, for which so far 
no other tools exist. 
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7. SUMMARY 



We summarize here what we take to be the main conclusions reached during the past 
two decades in the study of colour deconfinement and its manifestation in high energy 
nuclear collisions. 

• QCD, in particular through finite temperature lattice studies, predicts the onset 
of colour deconfinement (quark-gluon plasma formation) for strongly interacting 
matter of vanishing baryon number, once energy densities in the range 1-3 GeV 
are reached. 

• High energy nuclear collisions produce in their early stages sufficiently high energy 
densities for deconfinement; however, the produced systems expand and cool rapidly. 
Hence the crucial problem is the specification of probes which can unambiguously 
test the confinement/deconfinement status of the produced early medium. 

• For a direct test, the high density, short distance nature of the quark-gluon plasma 
requires hard probes which can distinguish between confined and deconfined media. 
Two such probes are: 

- quarkonium dissociation, which for fully formed physical resonances can only occur 
in a deconfined medium, and 

- jet quenching, manifested by a much higher energy loss of hard transverse jets in 
a quark-gluon plasma than in hadronic matter. 

• Charmonium suppression provides a particularly well-studied example, for which 
today there exist both an established theoretical framework and quite conclusive 
experimental results from nuclear collisions. 

- The hadroproduction of charmonium is understood as perturbative cc formation by 
interaction of partons whose distribution functions are determined in deep inelastic 
scattering; the colour neutralization of the cc pair can be described quantitatively 
in terms of colour evaporation. 

- Charmonium production in a confined medium is studied in p — A collisions, and 
the resulting A-dcpcndcnt decrease in production rate is understood in terms of 
pre-resonance absorption in standard nuclear matter. 

- Since about 60% of the observed J/ip^s are directly produced, while the rest 
comes from decay of higher cc excitations (mainly Xc), the onset of deconfinement 
in nucleus-nucleus collisions should lead to a suppression hierarchy. In a first step, 
essentially at the deconfinement point, the higher excitations arc suppressed. In a 
second step at clearly higher energy density, the smaller direct J/^'s are dissociated. 

- J/'i/' production in nucleus-nucleus collisions up to central S — U collisions shows 
only the normal pre-resonance absorption known from p—A collisions. While periph- 
eral Pb — Pb collisions also follow this pattern, there is with increasing centrality an 
abrupt first onset of some 30% further 'anomalous' suppression, followed by a later 
second further drop of some 20% in the observed J/ip suppression. This step-wise 
sudden onset of anomalous J/ ip suppression disagrees with all conventional hadronic 
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suppression models and is understandable only as the onset of colour deconfinement. 



The onset of deconfinement is supported by multiplicity data of soft hadronic sec- 
ondaries, which can be used to estimate the initial energy densities of the collision. 
With values of up to 3.7 GeV/fm^ in the inner region of central Pb — Pb collisions, 
the deconfinement range of energy densities is indeed attained. 

The relative strangeness abundance could provide another check of the initial energy 
density. If the strangeness suppression observed in p — p collisions of average asso- 
ciated multiplicity is due to the relatively low energy density of the early medium, 
it can be expected to end when in A — A interactions the initial energy density 
becomes high enough. There are indications that this is happening. 



J/ijj production experiments at the CERN-SPS [|107| , |108|| thus provide first evidence for 
the predicted new state of deconfined quarks and gluons. Future experiments at the RHIC 
in Brookhaven and the LHC at CERN, making use of charmonium, bottonium and jet 
production, are expected to extend these results and study the quark-gluon plasma in more 
detail and up to much higher energy densities. Soft probes, in particular the strangeness 
abundance, can serve to determine these initial energy densities quantitatively and study 
their effects on the final hadronic state of matter. 
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